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ABSTRACT 


A  coaxial  electromagnetic  shock  tube  has  been  used  to 
study  ionizing  shocks  propagating  through  hydrogen.  An  experi¬ 
ment  is  described  which  produced  plane  switch-cn  ionizing 
shocks.  Studies  were  conducted  covering  the  sub-Alfvenic, 
trans-Alfvenic  and  super -Alfvenic  regimes.  Data  is  presented 
that  verifies  the  existence  of  switch-on  fronts  in  all  these 
regimes,  and  demonstrates  the  significant  effect  of  the  ini¬ 
tial  electric  field  on  the  shock  jump  conditions.  The  results 
substantiate  the  theory  of  normal  ionizing  shocks  developed 
by  Taussig. 

A  simple  physical  criterion,  modeled  on  the  Chapman- 
Jouguet  (C-j)  hypothesis,  is  employed  to  choose  from  among 
the  several  theoretically  predicted  ionizing  waves  propagating 
at  a  given  speed.  This  criterion  has  strong  intuitive  support 
and  results  in  a  consistent  picture  of  ionizing  shock  propa¬ 
gation  that  is  in  general  agreement  with  our  experimental  ob¬ 
servations.  The  existence  of  a  small  homogeneous  plasma  sam¬ 
ple,  behind  the  ionising  shock  wave  and  before  the  expansion 
fan,  is  demonstrated.  Previous  failures  to  find  large  homo¬ 
geneous  plasma  samples  are  shown  to  be  consistent  with  the 
theory  of  ionizing  shock  waves  and  the  C-J  hypothesis. 
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The  present  theory  of  shock  phenomena  in  fluids  has  its 
roots  in  the  middle  of  the  nineteenth  century.  At  that  time, 
Riemann  and  Sarnshavr^^^  studied  long  wavelength,  large  an^)!!- 
tude  disturbances  propagating  in  media  described  by  non-line¬ 
ar  partial  differential  equations.  Their  studies  lead  to  the 
prediction  that  under  certain  circumstances  very  narrow  regions 
would  form  and  propagate  through  a  fluid.  Across  those  regions 
the  macroscopic  properties  of  the  fluid  were  expected  to  ex¬ 
hibit  large  scale  changes.  Those  narrow  spatial  regions  are 
what  we  know  today  as  shock  waves. 

At  the  turn  of  the  century,  Vielle^  conducted  the  first 
esqjerimental  study  of  shocks,  but  after  his  work  interest  un- 
e>?)lainably  declined.  Alnost  no  studies  were  conducted  during 
the  next  forty  years  with  the  exception  of  those  conducted  by 
Payman  and  his  associates^  in  the  “thirties."  The  hiatus  was 
abruptly  terminated  by  the  onset  of  the  Second  World  War,  The 
Allies  early  recognized  hov/  little  was  known  about  the  problems 
of  detonation,  explosion,  trans-sonic  flow,  and  super-sonic 
flow,  areas  that  are  all  intimately  connected  with  the  forma¬ 
tion  and  propagation  of  shock  fronts. 

During  the  war,  Bleakney  and  his  co-workers  at  Princeton 
started  research  into  those  areas.  As  a  result  of  their  stud¬ 
ies  the  diaphragm,  gas  driven,  shock  tube  became  the  basic  tool 
of  the  scientist  working  in  high  tenperature  gas  dynamics.  By 
the  early  1950* s,  Laporte  at  Michigan  and  Kantrowitz  and  his 
co-workers  at  Cornell  had  so  refined  the  shock  tube  that  they 
reached  the  practical  limit  of  its  ability  to  produce  high 
temperature  gases.  They  concluded  that  the  maximum  tenpera- 
ture  obtainable  in  the  conventional  gas  driven  shock  tube  would 
be  about  2  e.v. 
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But,  by  this  time,  the  controlled  thers»nuclear  program 
was  ui^er  way  and  it  called  for  laboratory  analysis  of  dense 
gases  at  hundreds,  and  even  thousands,  of  electron  volts.  By 
1950  de  Hoffman  and  Teller^  had  develop«3  a  theory  of  shtxrks 
in  such  immensely  hot  fluids.  It  was  based  on  the  laodel  of 
an  infinitely  conducting  fluid  immersed  in  nagnetic  and  electric 
fields  and  it  served  as  the  starting  point  for  many  subsequent 
theoretical  studies. 

In  the  early  1950* s,  Fowler  et.  al.  developed  a  magneti¬ 
cally  driven  shock  tube  t'lat  was  to  become  the  forerunner  of 
the  coaxial  electromagnetic  shock  tube. ®  Concurrently,  the 
theory  of  de  Hoffmcin  and  Teller  was  being  e3q>anded,  dissected 
and  reassembled.  By  1959  the  admissible  solutions  to  the 
magneto-hydrodynamic  (mhd)  shock  jump  equations  were  well  un¬ 
derstood  and  thoroughly  explored.  Questions  of  mhd  shock  wave 
stability  ware  being  discussed  and  advances  in  understanding 
were  being  made.  Experimental  studies,  however,  were  rela¬ 
tively  crude  and  inconclusive.  Most  experiments  involved 
strong  shock  waves  propagating  into  a  cold  gas  and  hence  were 
far  from  the  mhd  shock  wave  analyses  based  on  an  infinitely 
conducting  fluid.  In  1959^  an  attempt  was  made  to  generalize 
the  original  mhd  shock  theory  so  as  to  include  the  possibility 
of  a  jump  in  conductivity  across  the  shock  front.  This  approach 
was  taken  by  two  Russians^^  ®  who  investigated  the  mathemati¬ 
cal  properties  of  a  shock  propagating  into  a  gas  whose  conduc¬ 
tivity  was  assumed  zero  and  of  such  strength  that  the  gas  behind 
could  be  assumed  infinitely  conducting.  This  model  allows  the 
upstream  gas  to  support  an  arbitrary  electric  field  and  gives 
rise  to  the  type  of  discontinuity  that  will  occupy  us  in  this 
study  -  the  ionizing  shock.  Other  pioneering  work  in  this 
area  was  done  by  Gross  and  Kankel^^  and  Helliwell  and  Pack"“, 


in  19^.  Ifore  recent,  studies  will  be  described  in  the  body 
of  tSiis  renort.  With  this  introduction  we  turn  our  attention 
to  the  report  itself. 

Chapter  2  is  divide  into  tv?o  parts.  The  first  is  devo¬ 
ted  to  a  concise  e^roosition  of  the  ‘’classical”  theory  of  nsag— 
aeto-hydr odynami c  shocks,  i.e.  shocks  allowed  in  a  gas  wbose 
coiMSuccivity  is  infinite,  while  the  second  presents  ai»  e3q>lains 
the  most  inun^rtant  features  of  the  latest  theory  of  ionizing 
shocks. 

Chapter  3  bas  five  sections.  The  first  two  constitute 
a  conplete  report  on  the  construction,  operation  ai»  instru¬ 
mentation  of  a  coaxial  electromagnetic  shock  tube.  The  third 
eiqjlains  reduction  of  the  raw  data  to  a  form  convenient  for 
comparison  with  theory.  The  fourth  presents  a  proposal  that 
describes,  sinroly  ai^  consistently,  the  behavior  of  ionizing 
shocks.  The  proposal  is  based  on  the  theory  developed  in  the 
second  section  of  Chap,  2  and  on  the  data  in  the  third  section 
of  the  third  chapter.  The  reduced  data  is  compaured  with  the 
proposal  and  conclusions  are  drawn  in  the  fifth  section. 

Four  appendices  are  included  to  present  material  men¬ 
tioned,  but  not  treated,  in  the  text.  The  first  two  are  sum¬ 
maries  of  the  theory  of  shaped-pulse  capacitor  bamks  and  cur¬ 
rent  monitor  probes.  They  have  been  included  for  completeness  ai^ 
to  simplify  the  reader’s  approach  to  each  aurea.  The  third 
(Appendix  C)  develops  a  complete  picture  of  the  application 
of  coil  probes  to  the  measurement  of  magnetic  fields  in  moving 
media.  The  treatment  is  significant  in  that  it  points  up  the 
danger  of  using  an  oversimplified  view  of  such  probes  «md  fur¬ 
ther,  a  general  result  is  derived  that  is  applicable  to  all 
such  instrumentation.  Appendix  D  has  been  indluded  to  present 
additional  data  obtained  during  these  experiments.  Some  data 
on  ionizing  shock  formation  time  emd  ionizing  shock  thickness 
is  contained  therein. 


-3- 


I 


2.  HISTORY  AND  THEORY  OF  MAGNETO -HYDRODYNAMIC  SHOCKS 


ir.  2.1  "Classical"  Magneto-hydrodynamic  Shocks 

S'' 

:#  2.1.1  Theory 


it  is  generally  acknowledged  that  the  first  treatment 
of  magneto-hydrodynamic  (mhd)  shocks  was  the  work  of  de  Hoff¬ 
man  and  Teller^  in  1950»  However,  their  treatment  was  by 
no  means  exhaustive,  and  it  was  not  until  the  work  of  Bazer 
and  Ericson  in  1959  that  a  careful  study  was  done  of  the 
admissible  solutions  to  the  shock  jump  equations.  The  work 
of  Bazer  and  Ericson  leaned  most  heavily  on  two  earlier  papers. 

The  first  was  an  incisive  study  by  Friedrichs  ^  of  non-linear 
wave  motion.  He  was  the  first  to  recognize  the  "switch-on''  be¬ 
havior  of  one  class  of  shock  jump  solutions.  And  the  second,  which 
lead  directly  into  the  work  of  Bazer  and  Ericson,  was  a  study 
by  Friedrichs  and  Kranzer which  amplified  and  extended 
Friedrich's  earlier  analysis. 

At  the  same  time  that  these  papers  were  being  authored, 
independent  work  was  progressing  in  the  same  area.  Between 
1950  and  1959  important  work  was  completed  by  Lust^^  who  de¬ 
rived  directly  the  non-relativi Stic  mhd  shock  jump  equations 
(in  contrast  to  the  work  of  de  Hoffman  and  Teller  which  was 
done  within  a  relativistic  framework)  and  by  Lundquist’®  who 
was  the  first  to  treat  the  non-linear,  hyperbolic  system  of 
partial  differential  equations  that  govern  non-dissipative 
magneto-hydrodynamic  flow.  Even  as  far  back  as  19^9  ^nd 
1950'^^'®  there  v/ere  studies  showing  the  remarkable  feature 
that  characterizes  mhd  flow,  the  existence  of  three  distinct. 
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directionally  dependent,  small-disturbance  propagation  speeds. 
However,  these  earlier  studies  evolved  in  connection  with  a 
linearized  approach  to  the  mhd  problem  and  so  could  not  treat 
shock  phenomena. 

To  set  the  framework  for  subsequent  discussions  a  nec¬ 
essarily  brief  discussion  of  "classical"  mhd  shock  theory  will 
be  presented.  Here  the  word  "classical,"  far  from  meaning 
old  or  well  understood,  is  simply  being  used  to  distinguish 
mhd  shock  studies  done  in  infinitely  conducting  media  from 
the  more  recent  work  done  on  ionizing  shocks,  where  the  pre¬ 
shocked  gas  is  considered  to  be  completely  non-conducting 
(a  =  O),  and  only  the  post-shock  gas  is  assumed  to  be  in¬ 
finitely  conducting. 

The  appearance  of  shocks  in  any  real  fluid  is  bound  up 
with  the  non-linoar  nature  of  the  partial  differential  equa¬ 
tions  that  govern  the  flow.  In  the  study  of  infinitely  con¬ 
ducting  media  with  all  dissipative  coefficients  taken  as  zero, 
the  equations  form  a  system  of  eight  scalar  partial  differen¬ 
tial  equations  coupling  Maxwell’s  laws  with  the  fluid  equa¬ 
tions  that  speci:Qr  conservation  of  mass,  momentum  and  energy. 

The  equations,  often  called  the  Luiidquist  equations,*®  are: 

M  -  V  X  ( v  X  H)  =  0  (1) 

at 

+  v*V)p  +  pV’^'v  =  0  (2) 

+  v.V)v  +  lvp(p,s)  -  LL(V  X  H)x  H  =  0  (3) 

dt  p  p 

(—  +  v.7)s  ==  0  (4) 

dt 

where  H  is  the  magnetic  field,  v  is  the  fluid  velocity 
with  respect  to  any  fixed  reference  frame,  p  is  the  fluid 
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density;  d  is  the  fluid  pressure  considered  as  a  function 
of  the  two  thermodynamic  variables  p  and  S  (entropy), 
and  S  is  the  fluid  entropy. 

Equation  (l)  is  derived  from  Maxwell's  equations  with 
displacement  currents  neglected  (a  (~)^  ^  ^  1  approximation) 
and  a  (electrical  conductivity)  assumed  infinite.  Equation 
(2)  is  the  well-known  statement  of  mass  conservation.  Equa¬ 
tion  (3)  is  a  statement  of  conservation  of  momentum  in  the 
fluid.  It  is  derived  by  assuming: 

a)  that  V  X  H  =  j  so  that  the  j  x  B  term  in  the 
Lorentz  force  can  be  written  as  p,(  \/  x  H)  x  H: 

b)  that  the  Debye  shielding  distance  is  so  small  in 
comparison  to  ciny  flow  length  that  the  body  force  due  to 
interaction  of  the  electric  field  E  with  any  net  charge 
density  can  be  neglected ; 

c)  that  gradients  are  sufficiently  gentle  so  that  the 
scalar  pressure  term  dominates  the  viscous  force  term,  and 

d)  that  gravitational  effects  are  neglected. 

Equation  (  4)  states  that  the  system  is  reversible  and 
adiabatic  -  entropy  is  conserved.  This  equation  is  derived 
from  the  general  energy  conservation  equation  by  neglecting 
all  dissipative  effects.  It  is  obvious  that  Eq.  (4)  simply 
means  that  p  and  p  are  connected  through  the  adiabatic 
equation  of  state  p  =  A(s)p'^. 

A  system  of  equations  such  as  just  described  would 
seem  at  first  sight  to  be  so  restrictive  as  to  be  inapplicable. 
Just  the  contrary  is  the  case,  for  one  need  only  assume  that 
the  particle  distribution  function  in  the  fluid  is  "local" 
Maxwellian  to  prove  that  it  is  precisely  these  equations 
that  will  govern  mhd  fluid  flow.^°  The  achievement  of  a  "lo¬ 
cal"  Maxwellian  distribution  is  obviously  effected  through 
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collisions  and  so  the  validity  of  the  Lundquist  system  depends 
on  how  fast  and  within  what  length  scale  an  arbitrary  distri¬ 
bution  is  "Maxwellized, ”  Petsheck  and  Kantrowitz  have 
discussed  reasonable  criteria  for  estimating  the  range  of 
temperature  and  electron  density  within  which  the  Lundquist 
equations  are  valid.  They  argue  that  the  characteristic 
length,  L,  of  the  flow  field  must  satisfy  two  conditions: 


L  » 


3  X  10 


N. 


(5a) 


L  »  2/T‘ 


(5b) 


where  L  is  in  centimeters. 


is  teir^erature  in  e.v.  and 


is  the  electron  number  density  in  particles/cc.  If  5a 
and  5b  are  satisfied,  the  Lundquist  equations  can  be  assumed 
to  govern  the  behavior  of  the  flow  field.  The  inequalities 
5(a)  and  5(b)  represent  a  dimensional  analysis  of  the  assump¬ 
tions  used  in  deriving  the  Lundquist  equations.  For  example, 
if  we  assume  a  >]>  0  we  are  actually  requiring  that  the  mag¬ 


netic  Reynolds  number. 


R  =  u.avL, 
m 


be  much  greater  than  1, 


where  L  is,  as  before,  the  characteristic  length  of  the 
flow  field  and  v  is  the  characteristic  fluid  velocity. 

Further,  the  isentropic  assumption  is  equivalent  to  saying 
that  during  the  characteristic  flow  time,  t,  the  amount  of 
heat  added  to  the  fluid  via  conduction  is  much  less  than  the 
thermal  energy  of  the  plasma.  Thus,  we  require  k  «  pCLT 

l2  V  F 

where  k  is  the  heat  conduction  co-efficient,  and  is  the 
heat  capacity  at  constant  pressure.  In  a  similar  fashion, 
we  require  L  )>)>  Debye  length  to  insure  electrical  neutrality. 
Then,  using  kinetic  theory  expressions  for  k,  and  a,  5(a) 
and  5(b)  can  be  constructed. 
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If  an  attempt  is  made  to  construct  a  solution  to  the 
system  (l)  through  (4),  subject  to  explicit  boundary  conditions 
and  using  any  sensible  power  series  expansion,  one  finds  that 
there  exist  special  surfaces,  called  characteristic  surfaces, 
on  which  one  cannot  specify  initial  data  and  still  construct 
a  solution,  and  across  which  there  can  exist  jump  discontinu¬ 
ities  in  the  derivatives  of  the  dependent  variables.  This  is 
a  direct  result  of  the  hyperbolic  nature  of  the  system  of  equa- 

tionst®  In  addition,  as  is  well-known  in  ordinary  gas  dyna¬ 
mics,  there  is  a  critical  time  (after  the  initiation  of  some 

compressive  fluid  disturbance)  beyond  which  a  unique  continu¬ 
ation  of  a  smooth  solution  to  a  flow  problem  is  impossible. 
After  this  time  the  situation  is  no  longer  described  by  the 
Lundquist  equations,  and  finite  discontinuities  or  "jumps" 
will  appear  in  the  dependent  variables.  This  is  the  mathe¬ 
matical  origin  of  shock  waves. 

The  shock  problem  is  immensely  simplified  as  a  result 
of  the  conservation  form  of  the  differential  equations  that 
govern  the  continuous  flov?,  i.e,  as  we  have  seen  in  the  Lund¬ 
quist  equations,  each  equation  consists  of  the  sum  of  deriva¬ 
tives  of  functions  involving  the  dependent  variables.  From 
this  fact  it  can  be  shown*®  that  across  a  shock  the  dependent 
variables  are  related  by  a  system  of  elgebraic  equations  di¬ 
rectly  derivable  from  the  original  partial  differential  equa¬ 
tions  through  any  one  of  several  equivalent  approaches.  One 
may  use  a  generalized  Rankine-Hugoniot  relation,^®  or  a  siirple 
substitution  technique  first  used  by  Friedrichs,  or,  finally, 
the  direct  integration  approach  which  treats  the  shock  as  an 
infinitesimal  discontinuity  and  requires  integration  from  “he 
constant  state  ahead  of  the  shock  to  the  constant  state  be¬ 
hind.  Each  method  gives  the  same  set  of  algebraic  jump 

equations  relating  flow  properties  on  one  side  of  the  shock 
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to  flo.v  properties  on  the  other.  If  the  set  of  jump  equations 
is  viewed  from  the  rest  frame  of  the  shock  and,  if  it  is 
assumed  that  the  flow  is  one-dimensional  (in  the  x-direction), 
then  all  the  dependent  variables  are  functions  only  of  x 
and  t,  and  the  jump  equations  become: 

H  =  constant  A 


H  u  -  H  w 

II 

» 

c 

-  H  w 

A2 

Z8  2  X  2 

Z1  1 

X  1 

P2U2 

=  Pi  «1  = 

:  m 

A3 

p  u^  +  p  +  22 

=  p  u  ^ 

+  p  +  — ^ 

A4 

*2  2  2  2 

1  1 

^  2 

p  U  W  -  M-H  H 

^222  X  Z2 

II 

-  UH  H 
^X  Zl 

A5 

7-1  p. 
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u  ^+w  ^ 


(H  u  -H  w  )  = 
'  Z2  Z2  2  X  2 


m(-X-!i 
7-1  p, 


P.  \ 


u  ^  +  w  ^ 


-)  a6 


+  (H^  u  -H^w  ) 


m[S2(p,p)  -  S3^{p,p)j  >  0 


where  the  numerical  subscript  1  (2)  denotes  the  region 

ahead  of  (behind)  the  wave,  the  flow  velocity  vector  v  has 
components  (u, 0,w)  and  the  magnetic  field  vector  H  has  com¬ 
ponents  (h  ,0,H  ).  The  geometry  is  shown  in  Fig.  1.  The 
fact  that  the  flow  vector  and  the  magnetic  field  ahead  of  and 
behind  the  wave  can  be  restricted  to  the  x-z  plane  without  any 
consequent  loss  of  generality  is  proved  in  Bazer  and  Ericsoni^ 
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To  establish  the  basic  situation  with  which  we  will  be 
concerned  in  all  subsequent  discussions^  let  us  now  require 
that  the  fluid  ahead  of  the  shock  be  at  rest  in  the  laboratory 
reference  frame,  and  even  more  important,  that  H  =  0.  It 
is  this  latter  restriction  which  gives  rise  to  the  designa¬ 
tion  "normal"  shock  since,  in  this  case,  the  shock  front  is 
normal  to  the  magnetic  field  ahead  of  it.  By  virtue  of  the 
fact  that  the  shock  frame  has  been  assumed  in  motion  along 
the  x-cixis  we  get 


0  -  = 

u. 

.V 

+  U  or  u,  =  -U  ,  and 
s  i  s^ 

(6) 

0  =  LP 
w  = 

1 

w 

1 

(7) 

where  the  superscript 

LF 

means  that  the  quantity  is  to  be 

considered  in  the  laboratory  frame  and  U  is  defined  to  be 
the  velocity  of  the  shock  with  respect  to  the  laboratory 
frame.  With  these  substitutions  A1-A6  become  A1'-A6': 
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Assuming  that  we  know  u,  ,  PijPi>^v  >  system 

A1*-a6*  constitutes  a  set  of  6  equations  in  the  seven  unknowns; 


u=.h_,h 


X2^ 


,Uo>W_,p, 


and  pg.  The  system,  therefore,  determines 
a  one-parameter  feunily  of  solutions  for  the  downstream  states 
in  terms  of  some  given  upstream  state.  We  note  that  AJ  has 
not  been  considered  since  it  only  serves  to  separate  the  en¬ 
tropy  increasing  solutions  from  the  entire  set  of  solutions. 

A  ccniplete  discussion  of  all  the  admissible  solutions  to  A1’-a6' 
(excluding  stability  considerations)  is  given  in  Bazer  and 
Ericson.^  A  review  of  some  major  features  of  the  solutions 
will  be  presented  below.  First,  however,  a  fevr  general  re¬ 
marks  about  mhd  shocks  are  in  order. 


In  striking  contrast  to  the  situation  in  ordinary  gas 
dynamics,  the  mhd  jump  equations  admit  two  distinct  kinds 
of  shocks,  each  of  which  propagates  with  a  distinct  velocity, 

A  plausibility  argument  that  illustrates  this  behavior  pro¬ 
ceeds  as  follows.  In  conventional  coirpressible  fluid  theory 
a  shock  forms  in  the  following  way.  Let  an  infinitesimal, 
compressive  pressure  pulse  be  driven  into  a  cold  gas.  The 
disturbance  propagates  through  the  cold  gas  at  the  sound 
speed,  a^^ .  The  gas  behind  the  pulse  is  at  a  slightly  higher 
temperature  by  virtue  of  the  pulse* s  conpressive  nature.  The 
sound  speed  in  the  ccn^iressed  region  is  now  a^  a^^  since 
the  square  of  the  sound  speed  depends  solely  and  linecurly  on 
the  temperature.  If  a  second  compressive  pressure  pulse  is  later 
initiated  at  the  same  boundary  as  the  first  it  propagates 
into  the  slighted  heated  gas  at  a^  y  and,  therefore 
catches  up  with  the  first  disturbance.  Similarly  a  third 
pulse  catches  the  first  two,  and  a  fourth,  the  preceding  three. 

In  this  way,  a  sharp  spatial  discontinuity  is  established  in 
the  gas  (i.e.,  a  shock  is  formed).  The  ordinary  gas  dynamic 
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equations  yield  only  one  small -disturbance  propagation  speed, 
the  sound  speed,  and  hence  only  one  ki.nd  of  shock  can  be  pro¬ 
duced.  The  rahd  equations,  however  yield  three  different,  small 
amplitude  propagation  speeds.  They  cure  generally  denoted  by 

C-  C.  .  .  =  b  and  C  ,  .  and  they  are  characterized 

fast^  intermediate  n  slow- 

by  the  inequalities: 


^slow  ^  ^  ^  ^fast 
^slow  ^  ^n  —  ^fasf 


(8) 

(9) 


In  (8)  and  (9)  a  is  the  ordinary  sound  speed  in  the  gas 

and  the  “intermediate"  vave  soeed,  b  ,  often  called  trans- 

^  \^'Un  IH  i 

verse  wave  speed,  is  just  the  ordinauy  Alfven  speed,  b  =  — 2J — L 


with  I'Hl  replaced  by  the  component  of 
gation  direction  n. 


H  in  the  propa 


_ 


It  is  apparent  that  in  the  case  in  which  the  magnetic 
field  is  parallel  to  the  propagation  direction,  b^  =  b.  In 
the  same  case,  it  can  be  shown  there  is  a  degeneracy  that 
requires  c^  to  coincide  with  either  a  or  b  depending 
on  which  of  these  is  larger^  The  ratio  a/b,  then,  becomes 
cm  important  parameter  in  mhd  shock  theory.  From  now  on  we 
will  restrict  ourselves  to  the  situation  where  a/b  ^  1, 
since  it  is  only  in  this  case  that  the  shock  solutions  admit 
the  important  "switch-on"  wavet^  We  will  have  more  to  say 
about  this  later. 


In  exact  analogy  with  ordinary  gas  dynamics,  a  pressure 
puli.  .  propagating  at  any  of  the  three  small-disturbance  speeds 
might  be  expected  to  form  a  shock  and  therefore  we  would  ex¬ 
pect  to  find  three  different  shocks  possible  in  any  mhd  fluid 
flow.  However,  an  "intermediate"  or  transverse  small-distur- 
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bance  will  not  exhibit  any  steepening  tendency.  An  infini¬ 
tesimal  transverse  pulse  leaves  both  density  and  magnetic 
field  behind  it  unchamged.  Hence,  the  propagation  speed  be¬ 
hind  the  wave  is  identical  to  that  in  front  and  subseqpaent 
pulses  do  not  "catch"  the  initial  disturbance.  Large  aa5)li- 
tude  disturbances  of  the  transverse  variety  caui  be  propagated 
without  any  change  in  shape.  Only  fast  euid  slow  pressure 
pulses  exhibit  the  kind  of  steepening  behavior  that  chauracter- 
izes  a  pressure  pulse  propagating  in  cm  ordinary  gas  dynamic 
situation.  Thus,  it  is  the  fast  and  slow  waves  that  give  rise 
to  the  two  shock  modes  spoken  of  above.  Known  as  the  fast 
eind  the  slow  shock,  these  modes  derive  their  names,  quite 
obviously,  from  the  fact  that  the  former  is  formed  by  pressure 
pulses  traveling  at  the  fastest  of  the  three  small  distur¬ 
bance  speeds  and  the  latter  is  formed  from  pressure  pulses 
traveling  at  the  slowest  of  the  three  small-disturbance  speeds. 

We  now  turn  our  attention  directly  to  A1*-a6’,  Following 
Taussig^^  we  introduce  the  following  dimensionless  v«u:iables: 


"i 


lU-  j 

ill  w 


a.  w. 

r  ~  X 

,  a^  =  - w  =  - 

b  .  ^  b  . 

^  Xi  Xi 


•K  -  g- 
Xi 


i  =  1, 2  where 


b  . 

Xi 


=  1-^  H  ,  b  .  =  f-ii-  H 
VPi  ^  VPi  ^ 


Z1 


and  a 


i  =./ 


1^1 


V  Pi 


Substituting  intoAl*-A6’  the  ^stem  is  reduced  to  dimension¬ 
less  form.  Solving  for  u,  (assuming  for  convenience,  that 
-y  =  ^)  we  derive  a  fourth  order  ec[uation  for 

-  a/)[u^^  -  (4S3^  -  =  0  .  (10) 
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Equation  (lO)  obviously  has  four  roots,  one  of  which  is 
double: 

u^^  =  1^  a  double  root;  (lla 

u^^  =  (lib 

a  ^  ~  3a  ^).  (11c 

These  eu:e  plotted  in  Pig.  2.  With  these  solutions  for 
u^  in  terms  of  the  parciroeter  u^  (we  have  already  assumed 
a,  fixed  and  ^  l)  v;e  can  go  back  toAl’-A6'and  solv^e  ^stem- 
atically  for  all  the  shock  properties  in  region  2.  This  has 
been  done  by  Bazer  amd  Ericson,*^ 

The  solution  u^^  “  represents  the  trivial  "null" 

solution  to  the  jump  equations,  that  is,  all  the  shock  proper¬ 
ties  remain  unchanged  across  the  shock.  Though  trivial,  the 
null  solution  is  by  no  meams  uniirportant.  In  Fig,  2  it  repre¬ 
sents  the  dividing  line  between  regions  of  admissible  amd  in¬ 
admissible  solutions.  All  solutions  in  the  region  u^  ^ 
are  not  physically  realizable  because  they  imply  a  negative 
entropy  change  across  the  shock.  This  can  be  seen  immediately 
from  the  fact  that  for  u^,  Pe  ^  Pi*  expansive 

shocks  are  entropy  decreasing  shocks,  amd  are  therefore  not: 
allowed.  (For  proof  of  this  see  Bazer  and  Ericsoif'^  where 
it  is  also  proved  that  only  conpressive  shocks,  y  p,  , 
have  )>  ) .  in  addition,  the  null  solution  gives  rise 

to  an  important  class  of  solutions  in  the  theory  of  ionizing 
shocks  (see  Sec.  2.2.1). 

The  solution  u^^  =1  is  represented  by  the  straight 
line  BC  in  Fig,  2.  Note  that  to  the  left  of B  there  can  be 
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no  solutions  because  Pg  ^  Pi •  line  BC  is  the  switch- 

on  solution  that  has  already  been  mantioned  several  times  but 
has  not  been  discussed.  The  switch-on  wave  has  three  salient 
features.  First,  if  u^^  =  1  is  substituted  in  A1*-a6‘  and 
if  we  solve  for  H  ,  we  find  that  along  BC,  H  r  0.  A 
direct  result  of  this  fact  is  that  a  current  flows  through  the 
shock  and,  using  v  x  1^  -  the  current  density  is  of  precisely 
the  magnitude  necessary  to  produce  the  H  predicted.  The  quali- 

Z2 

tative  behavior  of  along  BC  is  easy  to  predict.  is  ini¬ 

tially  0  at  B,  as  it  must  be  since  B  is  on  the  null 
solution  where  properties  are  unchanged  across  the  wave.  It 
rises  to  a  maximum  at  some  point  along  BC  and  falls  to  0 
at  point  C.  Beyond  C,  H  is  imaginary  and  therefore, 
no  solutions  exist.  Thus,  across  a  switch-on  shock  there  is 
a  sudden  appearance  of  a  tangential  magnetic  field  ai»3  an 
associated  tangential  fluid  velocity,  deterodned  by  A2* , 

A  magnetic,  as  well  as  a  shear  discontinuity,  has  been  s%«ltched- 
on  behind  the  shock.  It  is  to  be  noted  that  only  the  switch- 
on  solution  (SW)  exhibits  this  behavior.  Since  this  shock 
has  no  ordinary  gas  dynamic  analog  (a  shear  discontinuity 
cannot  be  supported  in  ordinaary  fluid  flow)  it  has  been  the 
subject  of  intensive  theoretical  and  experimental  investiga¬ 
tion. 

The  second  feature  of  switch-on  shocks  is  that  the  flow 
field  behind  the  shock  has  a  normal  component  relative  to  the 
shock  which  is  exactly  the  normal  Alfven  velocity,  b^,  and 
which  remains  unchanged  as  the  shock  ^eed  increases.  And 
third,  the  switch-on  shock  divides  the  u^  -  u^  diagram  into 
two  distinct  regions,  the  slow  shock  region  <  1,  and  the 
fast  shock  region  Ug  ^  1,  where  the  designations  ”slow" 
and  "fast”  conform  to  the  definitions  of  Bazcr  and  Ericson,^^ 
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The  question  of  switch-on  stability  has  not  yet  been 
answered  in  full.  It  has  been  shown  that  the  switch-on  shock 
has  a  unique  structure  “^but  beyond  that  very  little  conclu¬ 
sive  work  is  available. 

The  solution  ^  been  plotted  as 

the  curve  C'ADCEco.  it  has  three  distinct  regions:  0‘A,  ADC, 
and  CE  «.  O' A  is  inadmissible  from  entropy  considerations. 

ADC  is  classed  as  a  slow  gas  shock  solution  and  CE  a 
fast  gas  shock  solution.  The  adjective  ‘'gas"  is  used  to  in¬ 
dicate  that  the  solutions  behave  somewhat  like  classical  gas 
dynamic  shocks  in  that  no  tangential  magnetic  field  or  velocity 
is  switched  on.  Though  i)e-ther  of  the  admissible  solutions 
exhibit  switch-on  behavior,  nevertheless  they  have  come  in 

i  Q 

for  extensive  review,  particularly  DC.  Taniuti  and  Jeffrey 
report  that  DC  is  a  non-evolutionary  solution,  i.e.,  it  i s 
a  solution  of  nhe  jump  equations  but  it  is  not  a  solution 
that  is  continuously  dependent  on  infinitesimal  changes  in 
boundary  conditions.  Gardner  and  Kruskal^'^  conclude  that  it 
is  not  stable,  v;hile  Petcheck  and  Kantrowits  have  argued 
for  the  interpretation  that  DC  is  an  extraneous  solution. 

All  the  arguments  suffer  from  -ne  following  criticism. 

If  we  must  exclude  DC,  we  may  ask  what  happens  in  a  piston 
problem  to  slow  shocks  that  are  driven  at  speeds  just  beyond 
those  corresponding  to  the  point  D,  Must  they  jump  to  a 
switch-on  type  of  discontinuity,  even  though  the  pifiton  need 
not  be  constrained  to  rotate?  Or,  if  not,  what  does  happen? 

The  questions  are  unanswered.  The  regions  AD  and  CEoo  arp» 
well  undersroodc  The  former  is  evolutionary  but  has  r't  yet 
been  shown  to  be  stable,*  The  latter  has  been  demonstrated 
to  be  both  stable  and  evolutionary,^® 
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2.1.2  Exper i men t  s 


! 


As  early  as  195^  strong  mhd  shocks  were  being  generated 
in  the  laboratory,  though  they  were  not  recognized  nox  studied 
as  such  at  the  time.  The  devices  then  used  were  known  as 
T-tubes  for  they  consisted  of  a  glass  vessel  in  the  form  of 
a  T.  Embedded  in  the  head  of  the  T  were  tv;o  electrodes 
betv/een  which  an  intense  electrical  discharge  v;as  initiated. 

By  bringing  the  return  current  path  near  the  discharge  {called 
back- strapping)  the  currenr  was  driven  by  the  Lorenfcz  force 
into  the  stem  of  the  T  forming  a  strong  shock  wave  in  the 
process.  By  1957  Kolb,^^  using  the  T-tube,  had  driven  a  deu-- 
terium  plasma  at  approximately  15  cm/pis,  graphically  demon¬ 
strating  the  immense  advantage  the  electromagnetic  driver  en¬ 
joys  over  the  conventional  pressure  driver.  However,  most 
of  the  T-tube  studies  were  concerned  with  producing  a  hot 
plasma  for  use  as  a  spectroscopic  source  and  little  work  was 
done  in  veri^ing  mhd  shock  jun^  properties. 

The  first  attempt  to  verify,  in  any  sense,  the  mhd  jump 
equations  was  made  by  Patrick  in  .1959^^  using  a  magnetic  annu¬ 
lar  shock  tube  (MAST).  Srncc  the  operation  of  the  Columbia 
coaxial  shock  tube  described  in  chapter  3  a  close  parallel 
to  the  MAST,  vie  n^ed  only  mention  that  the  MAST  provided 
cleaner,  more  reproducible  laboxatory  conditions  foi*  mhd  shock 
formation  than  the  T-tube. 

In  any  laboratory  situation  it  is  crucial  that  there  be 
availaole  an  accurate  theoretical  analysis  of  uhe  operation 
of  laboratory  devices,  and  the  mhd  shock  problem  is  no  excep¬ 
tion.  Without,  a  theory  of  electronuagneric  sho>ck  tube  opera- 
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tion.  it  would  be  a  hopeless  task  to  try  to  verify  the  shock 
jump  equations.  For  this  reason  much  wotk  has  been  done  on 
solving  the  shock  formation  problem  in  the  context  of  a  realis¬ 
tic  laboratory  boundary  value  problem.  The  earliest  t:heoreti- 
cal  model,  advanced  by  Garwin  and  Rosenbluth  in  195^j 
the  ’’snov^low”  model.  Still  in  use  today  as  a  good  qualita¬ 
tive  description  of  the  ‘'magnetic"  piston  (see  Sec.  3.1.1  for 
more  discussion  of  this  term),  the  theory  described  the  dr:.ve 
current  sheet  as  a  "snowplow"  pushing  gas  ahead  of  ic  and  es¬ 
tablishing  a  strong  shock.  In  1958,  Wright  and  Black  ex¬ 
tended  this  mode.l  by  including  acceleration  of  the  "snoivplow" 
and  discussing  the  time  that  must  elapse  before  steady  condi- 
conditions  prevail  in  any  electromagnetically  driven  shock 
tube.  Kemp  and  Petscheck^®  further  refined  the  shock  tubt 
problem.  They  allowed  gas  to  flow  through  the  snowplov?  by 
considering  a  rarefaction  fan  that  carried  the  drive  currant, 
but,  more  i  mportan  t,  they  motivated  thei..  '‘^'ork  with 
reference  to  the  experimental  design  of  Patrick^®  so  that  a 
device  and  a  theory  were  simultaneously  available.  The  analy¬ 
sis  predicted  that  between  the  shock  and  the  rarefaction  fan 
one  should  find,  as  one  v/ould  hope,  a  region  of  uniform  hot 
plasma,  i.e.,  precisely  the  uniform  region  behind  the  shock 
that  we  have  been  calling  region  2.  Patrick's  initial  work^® 
showed  that  reasonable  agreement  betv/een  observed  shock  velo¬ 
cities  and  those  predicted  by  the  Kemp  and  Petschek  model 
could  be  obtained.  Later  work  showed  marginal  evidence  of 
the  formation  of  a  homogeneous  gas  sample  {  hgs)  behind  a  strong 
mhd  shock,  and  also  some  of  the  magnetic  compression  effects 
that  occur  in  the  case  of  non-normal  fields. 

The  most  recent  attempts  to  verify  the  jump  equations, 
with  only  normal  fields  being  imposed,  have  not  been  conclu¬ 
sive.  Reiser^'  shows  evidence  of  switch-oir  .behavior  but  his 
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results  contain  9reat  scatter  and  are  open  to  differing  inter¬ 
pretations.  Hsyv;ood^“  states  that  for  the  normal  field  ccn= 
figuration,  the  Kerap  and  Petschek  model  does  not  describe 
shock  tube  behavior  since  neither  hgs  nor  switch-on  behavior 
could  be  observed  or  produced  in  his  experiment,  in  sum, 
through  January  1965^  there  has  been  no  complete  experimental 
study  of  the  ’'classical”  mhd  shock  problem,  let  alone  the 
special  normal  case. 


The  most  important  reason  for  the  lack  of  any  definitive 
experimental  studies  is  that,  the  classical  theory  requires 
that  any  mhd  shock  be  formed  in  an  already  infinit.ely  conduct¬ 
ing  fluid.  Since  such  a  fluid  is  not  readily  available,  the 
question  of  producing  one  becomes  pertinent.  One  way  tc  do 
so  would  be  to  take  a  cold  non—condueting  gas  (o  =  C),  drive 
a  strong  shock  through  it  and  make  it  infinitely  conducting. 
This  approach  lead:-,  quir.e  naturally,  to  the  study  of  the 
ionxzing  shock.  If  we  restrict  ourselves,  as  above,  to  only 
normal  magnetic  fields  we  hr.ve  normal  ionizing  shocks.  Their 
properties  and  their  theory  v;ill  now  he  presented. 
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2.2  Normal  Ionizing  Shocks 

2.2.1  Theory 

The  extension  of  "classical"  mhd  shock  theory  to  the 
case  where  the  fluid  is  assumed  to  have  finite  conductivity 
is  accomplished  by  the  trivial  addition  of  a  Poynting  v'ector 
term  to  the  energy  jump  equation. The  consequences  of  the 
addition  of  this  term,  however,  are  anything  but  trivial. 

Not  only  does  it  produce  a  whole  new  class  of  possible  solu¬ 
tions  to  the  jump  equations,  but  it  also  predicts  phenomena 
heretofore  totally  unexpected  and  unappreciated. 

By  a  normal  ionizing  shock  we  mean  that  we  restrict  our¬ 
selves  to  the  case  of  a  normal  magnetic  field,  a  stationary 
fluid  ahead  of  the  shock  in  the  laboratory  frame,  an  electric 
field  only  in  the  plane  of  the  shock  (E  ^0  everywhere) 
and  a  theoretical  model  in  which  the  gas  upstream  of  the  shock 
has  a  =  0  and  downstream  *3  ~  •  It  then  follov;s  quite 

simply,  that  the  jump  equations  governing  normal  ionizing 
shock  areAl'- a6’ with  the  addition  of  a  term  to  the 

left  hand  side  of  a6*  .  For  convenierce  we  will  -ewrite  these 
equations  as  B1-b6  ,  v;ith  all  quantities  in  the  shock  frame. 


H  =  constant 

(Bi) 

H  u_-K  w 
za  2  X  2 

=  0 

(B2) 

=  m 

(B3) 

^  2  .4.  „  -u 

(b4) 
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(B5) 


P  H 

r 2  Z  P  ^  X  22 


2  I  2 

•V  ’-^2  ''2 

Tz  ^  ^  -  Hx'a)  ^  V 

m(— i - +  -^) 

7-1  p  2 


22 


(b6) 


Using  B2  and  b6  reduces  this  system  to  Bl'-B5',  the 
system  that  will  be  used  from  now  on. 


H  =  constant 

X 

(Bi«) 

p  u  =  p.  a. 

(B2’) 

pu^4*p  +  -  PU^  +  P 

^  ^  ^ 

(B3‘) 

p  u  w  -  p.H  H  =0 

(B4’) 

u  “  +  w  ^  p.  a  ^ 

- 2-)  -f  Eya  =  m(-^  +  -^) 

2  -^7.^  7-1  p  2 

(B5’) 

It  is  to  be  noted  that  appears  without  subscript 

1  or  2.  This  is  a  direct  result  of  the  fact  that  Ey  is 
constant  across  the  shock.  Also  implicit  in  the  system  Bl’- 
B3’  is  the  fact  that  E.  has  been  assumed  to  be  of  the  form 
:C,  Ey,  0, ,  an  assumption  that  can  be  made  with  no  loss  of 
generality  for,  by  rotating  one’s  coordinate  system,  the 
vanishing  of  E  can  always  be  assured  (e  =  0  everywhere 
has  already  been  required). 

The  major  departure,  in  physical  terms,  of  the  ionizing 
shock  model  from  the  classical  model  is  the  ability  of  rhe 
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Upstteaiu  nOn-cOnduCting  gas  to  support  a  nonzero  eiectrio 
field  that  is  independent  of  flow  velocity  and  magnetic  field 
in  that  region.  Only  in  the  region  behind  the  shock  where 
a  =  00  must  we  require  that  E  =  -v  x  B.  Since  the  upstream 
electric  field  is  allowed  to  take  on  arbitrary  values,  we  ex¬ 
pect  to  appear  as  a  parameter  in  any  ionizing  shock 

theory.  This  is  precisely  the  case. 

The  particular  model  chosen  to  describe  ionizing  shocks 
was  adopted  because  it  permits  the  simplest  theoretical  treat¬ 
ment.  By  a  straightforward  exercise  in  the  use  of  the  Gali- 
lean  transformation  E‘  =  E  -f  v  x  B  one  can  easily  rewxite 
Ey  in  the  shock  frame  in  terms  of  the  downstream  variables 
(or,  since  Ey  is  a  parameter  it  might  be  more  appropriate 
to  think  in  terms  of  relating  the  downstream  variables  to 
Ey) .  The  relationship  is 


El-  =  |i{u  H 

r  2  Z2 


). 

a  x' 


(12) 


If  (12)  is  used  in  B5*  and  if  the  dimensionless  variables  al¬ 
ready  introduced  in  2.1,1  are  supplemented  by  an  additional 
dimensionless  variable  E  =  - = - .  then  E1’-3S’  can  be 

reduced  exactly  as  A1»-A6’  was  reduced  in  Sec.  2.1.1.  The 
generalization  of  (10)  becomes 


As  before  we  have  u  as  a 

~  2 

we  now  have  E  playing  the 
2.1  we  fix,  a  priori.  —  ^ 


function  of  u,  but  in  addition 
role  of  a  parameter.  (As  in  Sec. 
1.) 
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Equation  (13)  was  first  derived  and  discussed  in  full 
detail  by  Taussig,  although  the  actual  problem  was  first 
appreciated  by  Liubimov  and  Kulikovsky  in  1959.  ®  The 

substance  of  Taussig’s  work  will  be  presented  here.  There 
are  several  important  features  that  distinguish  normal  ioni¬ 
zing  shock  solutions  from  classical  mhd  shock  solutions,  ifost 
can  be  described  by  referring  to  Fig.  reprinted,  with 
additions,  from.  Taussig’ s  report.  The  heavy  black  lines  repre- 
sent  the  E  =  0  solutions  discussed  in  Sec.  2,1.1  and  pre¬ 
sented  in  Fig.  2.  For  a  given  E  )>  0  the  solutions  tc  (13) 
are  shown  as  the  light  curves  in  Fig.  3»  Three  features  are 
immediately  apparent.  First,  the  E  0  solutions  are  bounded 
by  the  E  =  0  solutions  and  this  is,  in  fact,  true  for  any 
E  )>  0.  Second,  if  E  is  specified,  there  are  ranges  for 
which  no  shock  solutions  atre  available,  e.g.,  the  region 
q,  <(  u^  \  q2  contains  no  real  solutions  to  Eq.  (13);  and 
also,  the  contours  that  describe  the  E  ^  0  solutions  are 
separate  and  distinct  as  is  evident  from  the  figure.  Third, 
contours  in  ABC  are  closed  curves  while  in  O’AO  and  E’BCE 
they  are  open  curves. 

Not  so  evident  is  the  fact,  proved  by  Taussig,  that  all 
solutions  for  E  ^  0  aire  nested  within  the  S  =  0  boundaury 
in  a  way  such  that,  for  any  two  values  of  E,  denoted  E^^ 
and  E^,  if  )>  E,  ,  then  the  E.,  contour  must  lie  within 
the  E.  contour.  Farther  arguments  can  be  advamced  to  prove 
that  not  only  aire  the  E  >  0  oolutions  nested,  but  there 

mm 

exists  a  critical  value  of  E,  denoted  by  E*,  such  that 
beyond  E*  there  exists  no  solutions  to  (13)  within  region 

mm 

ABC.  This  means  that  £is  we  increase  E  E*  the  closed  con¬ 
tours  in  ABC  shrink  to  a  point.  In  the  same  limit  of 
£  ^  E*,  ®  and  q,  0. 
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Each  contour  for  a  particular  E  )>  0  solution  is  divided 
into  several  sections  which  join  continuously.  The  sections 
are  labeled  K,  G,  and  SW  corresponding  to  the  portion  con¬ 
sidered  being  closest  to  the  null  (n).  gas  (g)  or  switch-on 
(SW)  solution  of  the  E  =  0  case,  respectively.  This  divi¬ 
sion  is  far  from  arbitrary.  It  corresponds  to  distinct  dif¬ 
ferences  in  the  jump  properties  across  the  shock.  One  deter¬ 
mines  the  jump  properties  as  before  by  going  back  to  the  full 
set  of  equations  Bl‘-B5*  and  systematically  solving  for  all 
variables  in  region  2,  given  u.  and  u„.  Two  important  re¬ 
sults  of  such  a  program  must  be  mentioned: 

1)  unambiguous,  switch-on  behavior  is  extended  to  all 
the  E  7  0  shock  solutions,  ana 

2)  the  null  branch  of  any  contour  becomes  a  full-fledged 
shock  in  contrast  to  the  "null  shock’’  in  the  E  =  0 
theory,  which  is  simply  a  mathematical  result  with 

no  physical  signi ficance. 

In  regard  to  the  first  of  these  points  v*e  note  that  switch- 
on  behavior  is  not  only  extended  to  slow  and  fast  shocks  but, 
even  mere  remarkably,  to  the  expansive  shocks  in  region  O’AO. 

This  brings  as  to  another  of  the  unejroected  features  of 
normal  ionizing  shocks  -  the  admissibility  of  expansive  shocks 
from  an  entropy  viewpoint. 

Consideration  of  the  Hugoniot  equation  for  normal  ioni¬ 
zing  shocks  leads  to  the  conclusion  that  in  ABC  all  solutions 
are  entropy  increasing,  while  in  O'AO  and  E’BCE  only  the 
solatitJHS  from  the  ’’bracket”  points,  ,  in  the  direction  of  the 
arrow  are  entropy  increasing  (see  Fig.  3)-  Thus,  as  stated  aixjve, 
at  least  a  sub  class  of  expansive  shoc3cs  is  admissible  as 
a  shock  solution  while,  as  already  noted  in  the  classical 
theory  ail  expansive  shocks  are  excluded. 
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The  problem  of  stability  and  evolutionarity  has  not  yet 
been  successfully  attacked  for  the  case  of  normal  ionizing 
shocks.  Only  transverse  ionizing  shocks  (i.e.,  ionizing 
shocks  where  H  is  parallel  to  the  shock  front)  have  been 
shown  to  be  evolutionary.  This  was  done  by  C.K.  Chu  in  1964.^'^ 
Chu’s  discussion  leads  to  a  remarkable  prediction.  In  order 
for  there  to  be  a  unique  structure  to  the  transverse  ionizing 
shock  the  electric  field  upstream  of  the  shock  cannot  be  arbi¬ 
trary,  as  has  been  assumed  throughout  this  discussion,  but 
must  be  altered  to  some  predictable  value  by  the  emission  of 
an  electromagnetic  wave  from  the  shock  front.  It  has  not  been 
established  that  this  result  is  carried  over  intact  to  tha 
case  of  the  normal  ionizing  shock  but  there  is  every  indica¬ 
tion  that  it  is.  Verification  of  this  point  awedts  further 
experimental  studies. 

Though  not  complete,  an  interesting  discussion  can  be 
advanced  to  argue  for  the  inadmissibility  of  all  E  0  so¬ 
lutions  in  the  region  Ug^  1,  for  the  case  of  the  normal  io¬ 
nizing  shock,  Recurning  to  34‘,  we  solve  for  v?^  and  sub¬ 
stitute  in  (12)  to  get; 


-  b^  )= 


liH  b^ 
Z2 


V. 


{^3^  -  1). 


(14) 


For  convenience,  let  us  assume  the  geometry  of  Fig.  1.  (The 
argument  is  identical  for  any  coordinate  ^stea. )  Now  ii 
and  b^_  are  always  0;  )>  1  by  hypothesis;  and  u^  )»  0 

for  the  geometry  chosen,  all  of  which  inplies  that  eii^  H 
are  always  of  the  same  sign.  But  this  would  seem  inpossible 
for  the  following  reason.  Since  is  continuous  across  the 
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shock  we  expect,  if  we  believe  j  =  oE,  to  have  j  parallel 
to  E.  But  this  i  will  give  rise  to  an  H  in  the  -z 
direction  contrary  to  the  statement  H  muse  be  greater 
than  0.  Taussig  includes  in  his  report  a  proof  by  C.K.  Chu 
that  attempts  to  substantiate  that  these  shocks  have  no  struc¬ 
ture  and  therefore  probably  do  not  exist.  From  such  arguments 
we  are  led  to  the  possibility  that  only  E  =  0  solutions  propa¬ 
gate  steadily  for  S  )>  1.  Since  the  method  by  which  ionizing 
shocks  are  initiated  would  seem  to  preclude  the  possibility  of 
having  E  =  0  ahead  of  the  wave  (see  Sec.  3-1)  the  question 
arises,  how  does  the  wave  adjust  to  the  E  =  0  requirement. 
Possibly  the  electromagnetic  wave  postulated  by  Chu  is  emitted 
from  the  shock  with  just  sufficient  strength  to  completely 
cancel  the  nonzero  electric  field  ahead  of  the  shock.  The 
ansv;er  to  this  question  can  only  be  provided  by  experiment. 

In  the  next  chapter  we  will  show  preliminary  evidence  that 
points  to  the  existence  of  E  /  0  solutions  in  apparent  dis¬ 
agreement  with  the  above  inadmissibility  aurgument. 

To  complete  the  theoretical  basis  for  an  experimental 
study  of  normal  ionizing  shocks  in  hydrogen,  Taussig  has  ex¬ 
tended  his  original  work  to  include  dissociation  and  ioniza¬ 
tion  across  the  shock  front. His  work  provides  a  con^juter 
program  whereby,  if  one  specifies  all  upstream  variables  in¬ 
cluding  shock  speed,  gas  corm>osition,  eleccric  field  and  nor¬ 
mal  magnetic  field,  then  the  downstream  state  is  determined, 
complete  with  percentage  dissociation,  percentage  ionization, 
temperature,  etc.  In  addition,  foilew’ing  Kenm*  and  Petschek 
exactly,  he  includes  integration  through  the  rarefaction  fan 
that  follows  the  shock.  Thus,  a  connjlete  theory  of  ionizing 
shock  properties  including  chemistry,  and  vcithin  the  context 
of  a  laboratory  situation,  is  available.  For  our  work,  such 
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a  refined  theory  is  not  necessary.  The  effects  of  chemistry 
serve  only  to  modify  slightly  the  non-chemistry  calculations. 

No  new  wave  solutions  result  from  inclusion  of  equilibrium 
chemistry,  and  hence  such  effects  will  be  iTC'Ortant  only  when 
the  downstream  variables  can  be  measured  to  within  a  few  per¬ 
cent.  However,  it  is  not  appropriate  to  neglect  the  fact  that 
if  one  deals  experimentally  with  hydrogen  (as  we  do),  ai^  if 
one  ionizes  the  hydrogen,  there  is  a  transition  from  diatomic 
to  monatomic  gas.  In  such  a  transition,  y  changes  from 
about  1.40  to  1.67.  This  effect  was  not  included  in  the  theory 
just  presented.  There,  y  was  assumed  ccnstauit.  If  one  does 
include  the  jump  in  y  in  the  jump  Eq.  B5',  the  form  aund 
all  the  features  of  the  solutions  remain  as  we  have  discussed 
them  -  only  the  numbers  change.  Ail  subsequent  work  in  this 
report  is  matched  to  the  theory  presented  in  this  chapter  with 
the  added  stipulation  that  there  be  a  jump  in  7,  across  tne 
shock,  from  1.40  to  1.67. 
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2,5,$  Experiments 


Frc'm  an  experimental  viewpoint,  ionisiing  shock  studies 
must  be  preliminary  to  any  studies  of  classical  mhd  shocks. 

This  is  a  point  that  has  been  too  often  overlooked.  It  is 
only  within  the  last  two  years  that  serious  work  has  begun  on 
producing  and  understanding  ionizing  shocks.  The  first  such 
studies  were  reported  in  I963  1964  by  Watson-Munro  et  al 

at  the  University  of  Sydney  in  Australia.  They  con¬ 

centrated  on  verifying  the  theoretical  results  of  Gross  and 
Kunkel,'’^  which  at  that  time  was  the  only  available  theoreti¬ 
cal  study  of  the  jump  properties  of  ionizing  fronts.  The 
Sydney  group  succeeded  in  verifying  the  wave  speed  predictions 
of  Kunkei  and  Gross. 

At  about  the  same  time  Vlases  reported  on  ionizing 
shocks  produced  in  an  inverse  pinch  e>rperiment. His 
work,  though  mostly  qualitative,  was  irnpcrtant  in  that  ha 
recognized  the  essential  differences  between  ionizing  waves, 
ordinary  gas  shocks  and  mhd  shocks.  In  mid~1965  Patrick  -...id 
Pugh,  reporting  on  their  work  with  ionizing  fronts,  report 
evidence  of  an  N  type  ionizing  front  in  a  coaxial  electro¬ 
magnetic  shock  tube.  However,  their  experiment  was  performed 
with  initial  magnetic  fields  in  the  plane  of  the  shock  front 
rather  than  normal  to  it. 

The  works  just  discussed  represent,  to  the  author’s 
>;nowledge,  all  the  experiments  conducted  specifically  to  study 
the  mechanics  of  mhd  ionizing  shock  fronts.  All  have  been 
mainly  concerned  v;ith  measuring  wave  speeds  and  mapping  out  de¬ 
vice  operation  and  none  has  been  significantly  extended  beyond 
that.  In  the  next  chapter  an  experiment  is  described  for  which 

t< 
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the  experimsnttil  e'/idlcr.cc  fits  a  simple  explanation  of  ioni" 
zing  shock  behavior.  The  proposal  that  is  advanced  and  veri¬ 
fied  is  inductively  constructed  from  the  theory  of  normal 
ionizing  shocks  just  presented,  from  experimental  evidence 
obtained,  and  from  analogy  with  classical  gas  dynamic  pheno¬ 
mena.  We  will  see  that  the  complete  range  of  shock  velocities 
can  be  mapped  by  a  single  curve  that  has  strong  intuitive 
support.  Further,  the  experiment  verifies  two  of  the  most 
unusual  aspects  of  the  ionizing  shock  front  -  the  appearance 
of  the  switch-on  shock  and  the  extension  of  the  switch-'on 
shock  to  the  slow  shock  regime. 


3.  AN  EXPERIMEiSrr  TO  VERIfY  TfiE  TIIKORY 
OP  NORT-IAL  IONIZING  SHOCKS 


3>1  Apparatus 

3.1.0  Introduction 


hn  experimental  verification  of  the  nor.naJ  ionizing 
shock  theory  presented  in  2.2.1  requires  a  device  capable  of 
generating  an  essentially  one-diraensional  shock,  propagating 
thr'ough  hydrogen,  whose  velocity  lies  in  the  range  10®  to 
greater  than  10^  cm/sec.  In  this  range  a  hydrogen  plasma 
can  be  produced  with  any  desired  degree  of  dissociation  and/ 
or  ionization. 

Attempts  to  attain  such  experimental  conditions  have 
been  made  with  conventional  pressure  driven  shock  tubes.  How¬ 
ever,  vath  these  devices  there  is  an  upper  lii.dt  of  about 
2  e.v.  for  the  temperature  that  can  be  produc'sd  in  the  shocked 
gas.  A  device  that,  seems  to  hold  great  promise  as  a  research 
tool  in  the  area  of  strong  ionizing  shocks  is  the  coaxial 
electromagnetic  shock  tube.  The  basic  principle  of  its  op¬ 
eration  is  simple.  One  discharges  a  large  amount  of  current 
from  one  electrode  to  an  adjacent  electrode,  and  the  seif- 
induced  magnetic  field  is  used  to  propel  this  current  into  a 
cold,  low-pressure  gas  region.  The  coaxial  geometry  requires 
that  the  discharge  be  across  the  annulus  between  cylinders 
and  so,  if  the  return  ground  is  carefully  arranged,  a  toroi¬ 
dal  current  flow  configuration  is  formed,  and  ail  the  self- 
induced  magnetic  field  is  confined  within  it.  The  Lorentz 
force  on  the  free- surface  of  the  current  flow  propels  that 
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surface  down  the  region  between  the  cylinders.  This  mechan¬ 
ism  is  analogous  to  forcing  a  piston  into  a  cylinder  and, 
like  its  mechanical  counterpart,  it  sets  up  a  shock  that  propa¬ 
gates  ahead  of  it.  Since  the  "piston”  can  be  driven  very 
fast,  the  shock  that  travels  ahead  of  it  can  generate  a  very 
high  temperature  plasma  ( see  Fig.  4  ) . 

The  experimental  apparatus  consists  ofs 

1. )  a  one  meter  long,  seven-centimeter  mean  diameter 

coaxial  shock  tube  with  a  one-centimeter  annulus; 

2. )  a  sclenoidai  magnet  capable  of  generating  longitu¬ 

dinal  d.c.  magnetic  fields  from  500  gauss  to  13.8 

kiiogaass; 

3. )  a  40, 000  joule  capacitor  bank  designed  to  deliver 

a  trapezoidal,  l4^s  pulse  of  150,000  amperes  into 

.  .2  ohm  load; 

4. )  an  air-gap  switch  for  switching  the  bank  into  the 

shock  tube; 

5. )  instrumentation  to  measure  shock  speeds,  induced 

magnetic  fields,  and  electric  fields; 

6.  )  control  circuitry  for  firing  the  shock  tube  and 

timing  the  diagnostic  devices. 


3,1,1  The  Shock  Tube 


The  shock  tube  is  basically  a  simple  device  consisting  of 
a  three  inch  i.d,  stainless  steel  cylinder,  one  meter  long. 

An  Inner  aluminum  cylinder  of  2^  inches  o.d,  is  mounted  co- 
cixial  with  the  stainless  outer  one.  One  end  of  the  inner 
cylinder  is  provided  with  a  copper  knife-edge  so  that  a  pre¬ 
ferred  breakdown  path  is  provided  when  the  capacitor  bank 


voltage  is  applied  between  the  inner  and  cuter  cylinders. 

The  knife-edge  is  mounted  on  a  |  inch  copper  rod  that  feeds 
through  a  piexi -glass  end  piece  to  one  end  of  the  air-gap 
switch  (see  Fig. 5  ).  The  air-gap  switch  then  connects  to 

the  high  voltage  side  of  che  capacitor  bank. 

The  knife-edge  is  positioned  on  the  inner  conductor  so 
that  after  assembly  inside  the  solenoidai  magnet  it  is  in  a  homo¬ 
genous  region  of  the  magnetic  field  (see  Fig. 6  ).  In  addi¬ 

tion,  to  provide  a  well-defined  boundary,  from  which  the 
shock  wave  is  launched,  a  flat  glass  plate  is  positioned  just 
behind  the  knife  edge  (see  Fig.  7  ).  It  was  found  necessary 
to  insulate  the  copper  feed,  and  the  last  9  inches  of  the 
inner  cylinder,  lest  the  discharge  be  initiated  there  instead 
of  at  the  knife  edge.  Those  areas  were  wound  with  thin  ny- 
lar  strips  and  teflon  tape,  epoxied  to  the  metal  to  insure 
a  homogeneous  seal. 

To  provide  an  experimental  situation  where  one  can  vai^ 
the  pressure  of  the  pre-shocked  gas  in  the  shock  tube  a  Cen- 
co  ti'pe  VMF  20,  20  1/ s  diffusion  pump,  backed  by  a  Welch  2 
CPM  fore-pump  is  attached  to  the  shock  tube.  This  system 
gives  a  base  pressure  of  approximately  lOp.  of  mercury  in  the 
working  region.  This  relatively  high  base  pressure  would 
seem  to  indicate  a  sizable  leak  in  the  si'stem,  but  leak  test¬ 
ing  with  a  mass-spectrometer  proved  otherwise.  The  glass 
backup  plate  behind  the  Icnife  edge  sc  throttles  pump-out  that 
although  a  base  pressure  of  5  x  10  ^  mm  of  mercury'^  is  obtained 
near  the  pump-out  port,  the  vacuum  system  cannot  reduce,  in 
any  reasonable  time,  the  base  pressure  in  the  working  area. 

The  base  pressure  near  the  pump-out  port  is  measured  by  a 
Veeco  Cold  Cathode  Discharge  Gauge  while  the  pressure  in  the 
working  region  is  monitored  by  a  Veeco  Thermocouple  Gauge 


attached  to  the  far  end  of  the  shock  tube.  The  thermocouple 
gauge  is  particularly  useful  because  the  system  is  generally 
filled  with  cold  hydrogen  to  about  200m.  of  mercury,  and  the 
thermocouple  gauge  can  measure  this  quite  easily. 

The  working  length  of  the  shock  tube  is  approximately 
75  centimeters  long,  and  consists  of  that  portion  of  the 
shock  tube  that  is  xn  a  relatively  homogeneous  region  of  the 
longitudinal  magnetic  field  ( see  Fig.  6  ) . 


3.1.2  The  Magnetic  Field 


The  magnetic  field  is  provided  by  36,  l.i”  wide  water- 
cooled,  copper  panccikes  potted  in  epoxy.  Designed  by  the 
Matterhorn  thermonuclear  group  at  Princeton  University,  these 
coils  are  structurally  able  to  withstand  peak  fields  of  50 
kilogauss.  They  produce  fields  homogeneous  to  +  over 

the  central  25  cm  inches,  and  ^  1^  over  the  central  33  eras, 

A  computer  program,  written  at  Columbia,  was  used  to  calculate 
the  field  produced.  Fig.  (6  )  shows  the  con^juter  calcula¬ 

tion  of  the  axial  and  radial  fields  as  a  function  of  magnet 
length.  Knife-edge  position  and  shock  tube  working  length 
are  also  indicated.  The  axial  field  was  measured  using  a 
hall  probe  gauss-meter  and  agreement  with  the  theoretical 
field  strength  was  at  least  as  good  as  the  gauss  meter  accu¬ 
racy  of  +  1%, 

The  magnet  is  powered  by  a  200  kilowatt,  three-phase, 
rectifier  with  less  than  1^  peak  to  peak  ripple.  The  recti¬ 
fier  was  designed  for  10-second  operation  at  300^  overload 
and  in  this  mode  of  operation  can  deliver  2200  aitperes. 

This  enables  the  magnet  to  provide  longitudinal  field  continu¬ 
ously  variable  from  about  500  gauss  to  almost  13.8  kilogauss. 
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3.1.3  The  Capacitor  Bank 


A  shaped-pulse,  40, 000  joule  capacitor  bank  provides  the 
energy  to  drive  the  ’’magnetic"  piston  in  the  shock  tube.  The 
bank  consists  of  a  total  of  18  capacitors  arranged  in  four 
parallel  impedance  meshes.  Each  mesh  is  equivalent  to  a  ca¬ 
pacitor  and  a  series  inductor.  This  particular  configuration, 
called  a  Guillemin  Type  C  network,  was  chosen  because  of  de¬ 
sign  simplicity,  ease  of  assembly  and  shaped-pulse  capability. 
Fig.  (8  )  is  a  schematic  of  the  bank,  and  gives  all  the  per¬ 
tinent  design  facts.  When  charged  to  50, 000  volts,  the  bank 
is  designed  to  deliver  a  1.0  M-s  rise-time  pulse  of  150,000 
amperes  into  a  matched  load  of  .2  ohms  for  approximately  l4 
microseconds.  However,  any  external  inductance  in  the  load 
badly  degrades  the  ideal  pulse  shape.  In  practice  we  get  a 
pulse  shape  more  nearly  sinusoidal.  Figure  (9  )  is  a  typi¬ 
cal  current  pulse  trace.  Appendix  (A)  has  been  included  to 
present  a  more  detailed  discussion  of  the  theory  behind  all 
shaped-pulse  capacitor  banks,  and  to  explain  the  motivation 
behind  the  Type  C  network. 

The  capacitor  bank  is  charged  to  its  operating  voltage 
by  a  Del  Electronics  (TC  series)  80kV,  ^xaa.  power  supply. 

This  power  supply  has  been  provided  with  automatic  shi  uotfs 
to  prevent  overcharging  the  ^OkV  capacitors,  and  equipped 
with  a  special  dump  tank  that  will  discharge  the  fully  charged 
bank  in  under  10  seconds. 

The  capacitor  bank  is  housed  inside  a  sturdy,  shielded, 
solid  wooden  enclosure  (Fig.  10}.  The  output  end  of  the  bank 
is  a  parallel  plate  transmission  line  whose  high-voltage  side 
is  connected  to  an  air  gap  switch. 
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3.1.4  Air-Gap  Switch 


The  air-gap  switch  separates  the  capacitor  bank  from 
the  copper  feed-through  attached  to  the  knife-edge  on  the 
inner  cylinder  of  the  shock  tube.  The  switch  is  most  easily 
described  by  referring  to  Fig.  (  11  ) .  The  solid  brass  anode 
attached  to  the  cheirged  capacitor  bank  is  a  flattened  hemi¬ 
sphere  adjusted  to  a  distance,  d,  from  the  stainless  steel 
cathode  that  is  sufficient  to  hold  off  the  capacitor  bank 
voltage.  It  is  important  to  note  that  some  d.c.  path  must 
be  established  from  cathode  to  ground  or  else  the  bank  vol¬ 
tage  divides  across  both  switch  and  shock  tube;  making  it 
difficult  to  fire  the  switch.  This  can  be  accomplished  ei¬ 
ther  by  connecting  a  large  (^10  Megohm)  resistor  between 
cathode  and  ground  or  by  pre-ionization  in  the  shock  tube. 

The  latter  method  is  the  one  presently  being  employed  because 
ic  serves  a  dual  purpose.  Discussion  of  this  is  deferred 
until  later. 

The  switch  is  fired  by  establishing  a  secondary  break¬ 
down  between  the  isolated  teaser  electrode  and  the  cathode. 

The  shock  formed  at  this  secondary  breakdown  carries  ions 
and  electrons  to  the  anode  and  establishes  a  low  impedance 
pack  between  anode  and  cathode.  The  capacitor  bank  voltage 
now  appears  at  the  knife-edge  on  the  inner  conductor  cuid  elec¬ 
trical  breakdown  occurs  there  as  desired. 

The  teaser  circuit  consists  of  a  10  kilovolt;  .125  M-f 
capacitor  charged  by  a  30  kilovolt,  5  nia  supply  through  a 
chcurging  resistor.  A  25  kilovolt  ignition  switch  holds  the 
voltage  until  firing  is  desired.  When  the  switch  is  to  be 
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firea,  a  I5OO  volt,  pulse  is  supplied  through  the  pulse  trans¬ 
former  between  ignitor  and  cathode  of  the  ignitrcn.  The  ig- 
nitron  breaks  down,  the  teaser  rises  to  10  kilovolts,  the 
teaser  breaks  down  to  the  cathode,  and  the  anode  then  breaks 
down  to  the  cathode.  Figure  (  12  )  shows  the  air-gap  switch 
in  normal  operating  position. 

Probably  the  most  inroortant  and  most  difficult  technical 
aspect  of  coaxial  shock  tube  technology  is  to  insure  that  af¬ 
ter  one  fires  the  capacitor  bank  the  dischcirge  at  the  knife- 
edge  on  the  inner  conductor  is  uniform  around  the  entire 
knife-edge.  Only  then  can  one  produce  the  "magnetic*'  piston 
geometry  which  is  the  key  tc  the  production  of  strong  ioni¬ 
sing  shocks. 

The  method  finally  decided  upon  in  this  experiment  was 
utilization  of  the  strong,  external,  longitudinal  magnetic 
field,  and  r.f.  pre-ionization  confined  to  the  immediate 
neighborhood  of  the  knife-edge. 


3.1.5  Pre-ionization:  Circuitry  and  Technique 


A  30  Mc/sec,  150  watt,  push-pull,  tuned  grid,  tuned-plate 
oscillator  was  designed  and  built  to  couple  symmetrically  to 
the  region  around  the  knife-edge.  Figure  (  13  )  gives  the 
oscillator  schematic  and  a  detail  of  the  coupling  arrangement. 
Localized  breakdown  was  initiated  in  hydrogen  at  a  pressure 
of  200  microns  of  mercury  by  imposing  a  longitudinal  magnetic 
field  of  greater  than  about  800  gauss,  and  then  pulsing  the 
oscillator  into  operation  at  about  1000  volts  and  100  ma. 
After  r.f.  breakdown  occurs  the  oscillator  level  is  reduced 
to  approximately  300  volts  and  10  to  15  ma.  The  application 
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of  the  longitudinal  magnetic  field  (i.e.,  the  crossed— field 
configuration)  vras  found  necessary  in  producing  uniform  r.f, 
breakdown.  Once  uniform  r.f.  breakdown  has  been  obtained, 
the  main  bank  discharge  can  be  initiated  and  it  breaks  down 
uniformly  around  the  knife-edge. 

Recalling  the  need  for  a  o.c.  path  from  air  gap  switch 
cath'^c  e  to  ground  (Sec.  it  is  obvious  that  the  pre¬ 

ionizer  performs  this  service  admirably. 


3.1.6  Control  Circuitry 


The  entire  experiment  is  controlled  from  a  single  panel 
where  one  push  button  can  fire  the  entire  apparatus  and  ini¬ 
tiate  ail  necessaury  timing  circuits.  The  basic  cinaing  cir¬ 
cuit  is  a  Pegram  Electronics  Laboratory  Dual-Delay  Module 
which  is,  sinroly,  a  variable  positive  or  negative  square 
pulse  generator.  Bach  unit  triggers  only  from  a  negative 
going  pulse,  and  this  feature  enables  one  to  set  up  aiy  num¬ 
ber  of  separate  pulses  of  either  polarity.  For  exanrole,  trig¬ 
ger  one  unit  and  obtain  a  positive  square  pulse  of  duration 
T  seconds.  Feed  this  into  the  trigger  of  the  second  unit, 
cuid  since  this  second  unj t  only  triggers  on  negative  going 
pulses  ir  will  not  trigger  until  the  trailing  edge  of  the 
first  pulse  arrives.  Thus  the  second  unit  dees  not  fire  un¬ 
til  T  seconds  after  the  first;  at  which  time  it  prouuces 
ai-  independent  pulse  of  either  poieurity  with  any  pre-set  du¬ 
ration.  Series  operation  in  this  manner  with  several  modules 
thus  allows  the  establishment  of  a  well-defined  time  base  with 
pulses  at  any  given  point  after  t  =  C,  The  t  =  C  point 
is  established  from  a  simple  blocking  oscillator  puls<»r  which 
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provides  a  negative  signal  sufficient  to  trigger  the  first 
dual-delay  module.  Associated  v;ith  the  timing  circuit  are 
two  high  voltage  pulsers.  One  is  a  fixed  duration,  transis¬ 
torized  pulser  capable  of  delivering  a  .2  microsecond,  600 
volt  pulse  into  30  feet  of  open  circuited  50  ohm  cable.  The 
ether  is  a  thyratron  circuit  which  discharges  a  1500  volt  capa¬ 
citor  into  the  1:1  transformer  in  the  air-gap  switch  circuit 
(Fig.  lo;*. 


3.1.7  Operation 


The  operation  of  the  shock  tube  consists  of  leaking  ul¬ 
tra-pure  hydrogen  into  the  annulus  through  a  \rariable  leak 
(see  Fig.  5  )  until  the  desired  operating  pressure  is  obtained 
as  monitored  on  the  thermocouple  gauge  (Fig.  5  )•  The  mag¬ 
netic  field  is  then  set  at  the  desired  value  and  pre-ioniza- 
tion  established  about  the  knife-edge.  It  has  been  observed 
that  for  reproduceable  operation  of  the  shock  tube  it  is 
necessary  that  once  pre-ionization  has  been  effected,  the 
power  level  into  the  r.f.  oscillator  be  reduced  to  a  minimum 
consistent  with  uniform  pre-ionization  (usually  this  is  of 
the  order  of  3~5  watts  plate  ..over,  whereas  breakdown  usually 
requires  about  150  watts) . 

The  capacitor  baiik  is  charged  to  the  voltage  that  will 
give  the  desired  drive  current  and  the  shock  tube  is  then 
ready  for  firing.  The  firing  sequence  is  initiated  by  a  push 
button  :n  the  blocking  oscillator  circuit  that  produces  a 
negative  pulse.  This  pulse  fires  the  first  dual-delay  module. 
The  output  from  this  module  then  triggers  all  the  oscillo¬ 
scope  time  bases  so  that  a  common  time  base  is  established. 
This  fii'st  pulse  also  fires  a  second  module  whose  output 
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triggers  the  6OO  volt,  fixed  duration,  high  voltage  pulser. 
Its  600  volt  pulse  is  delivered  to  the  the  thryatrcn  pulser, 
triggers  it,  and  the  thryatron  output  triggers  the  ignitron 
in  the  teaser  circuit  of  the  air-gap  switch.  The  air-gap 
switch  fires  and  the  bank  is  discharged  through  the  shock 
tube.  After  firing  the  ^stem  is  flushed  and  made  ready  for 
another  shot.  The  entire  timing  sequence  is  indicated  in  Fig. 
(14  ). 


. .  ' 


i.'d  Shock  Tube  Instrumentation 
3.2.0  Introduction 


Shock  tube  instrumentation  consists  of;  l)  a  simple 
magnetic  pickup  coil  for  measuring  total  drive  current  de¬ 
livered  by  the  capacitor  bank;  2)  photomultipliers  and  asso¬ 
ciated  circuitry  for  monitoring  the  velocity  of  the  luminous 
ionizing  front;  3)  magnetic  field  probes  for  detecting  the 
self-induced  field  in  and  behind  the  ionizing  front;  and 
4)  electric  field  probes  for  measuring  the  electric  field 
ahead  of  and  behind  the  shock. 


3.2.1  Current  Monitoring  Loop 


The  instantaneous  current  being  delivered  by  the  capa¬ 
citor  bank  is  determined  by  integrating  the  EMf-  (t)  observed 
between  the  terminals  of  a  coil  due  to  the  change  in  flux 
produced  by  the  time-varying  current. 

A  single,  square  loop,  carefully  shielded,  is  positioned 
along  the  m.ain  feed  into  the  shock  tube.  Its  terminals  are 
connected  to  a  50  ohm  cable  which  is  terminated  i  a  50  ohm, 
10:1  attenuator.  The  output  voltage  is  fed  into  a  Tektronix, 
Type  "o”  Operational  Amplifier  plug-in  unit  operating  as  an 
integrating  circuit.  The  output  of  the  operational  integrator 
is  l(t)  and  is  displayed  directly  on  the  scope.  Figs.  (15) 
and  ( l6)  show  the  loop,  its  construction,  its  position  on  the 
main  current  feed  and  the  circuit  used  to  measure  l(t).  The 
theory  of  the  probe  is  presented  in  Appendix  B. 
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3.2.2  Photomultipliers  and  Associated  Circuitry 


Photomultipliers  are  used  to  monitor  the  light  output 
from  the  shock  tube  at  several  different  ports  along  the 
shock  tube  length.  Light  is  delivered  to  a  photomultiplier 
through  a  fiber  optic  light  pipe.  The  light  pipe  is  connected 
to  a  reflector  mounted  above  a  l/l6'*  hole  drilled  through 
the  shock  tube  wall.  Port  details  are  shown  in  Fig.  (17b). 
Figure  (17a)  shov.’s  a  p’  tomultiplier  port  in  relation  to  the 
entire  shock  tube  assembly. 

RCA  93IA  photomultipliers  were  used  v;ith  standard  re¬ 
sistive  divider  networks  to  provide  voltage  to  the  10  stages. 
The  last  few  stages  are  provided  with  parallel  capacitors  to 
enhance  the  pulse  response  of  the  tubes. 

The  output  of  each  photomultiplier  is  connected  to  an 
emitter-follower  circuit  which  is  capable  of  driving  the  ter¬ 
minated  50  ohm  cable  that  brings  the  signal  to  an  oscillo¬ 
scope.  Figure  (  I8  )  is  the  complete  photomultiplier  circuit. 

3.2.3  Magnetic  Field  Probes 


The  theory  of  ionizing  shocks  distinguishes  three  regions 
of  interest  differentiated  by  the  value  of  the  Alfvenic  Mach 
Number,  M  In  the  sub-Alfvenic  region  (M  <(  l)  and  the 
trans-Alfvenic  region  (l\M  <(~2)  the  drive  current  di- 
vides  between  the  shock  front  and  the  current  carrying  rare¬ 
faction  fan.  In  the  super -Alfvenic  region  (M  )>  ~  2)  the 
drive  current  may  or  may  not  be  completely  confined  to  the 

rarefaction  fan,  depending  upon  whether  or  not  all  the  E  )>  0 

/ 

solutions  are  inadmissible.  Therefore,  the  sub-Alfvenic  and 
trans-Alfvenic  regions  v/ill  exhibit  sv/itch-on  behavior  while 
for  a  super -Alfvenic  wave  we  have  no  definite  theoretical  pre¬ 
diction. 
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It  is  apparent  that  measurement  of  the  seif-induced 
magnetic  field  holds  the  key  to  determination  of  switch-on 
like  ionizing  shocks.  In  the  switch-on  regime  the  magnetic 
field  profile  rises  through  the  shock,  remains  constant  in 
the  post-shock  gas,  and  then  rises  to  a  final  value  through 
the  rarefaction  fan.  In  the  super -Alfvenic  region,  the  mag¬ 
netic  field  profile  can  give  the  dramatic  evidence  of  separa¬ 
tion  of  shock  front  and  drive  current  that  is  necessary  to 
support  or  disprove  the  argument  about  inadmissibility. 

A  complete  theory  of  the  measurement  of  self-induced 
magnetic  fields  in  plasmas  will  be  deferred  to  the  next  chap¬ 
ter  and  to  an  appendix.  The  probe  is  designed  to  utilize 
Faraday's  law  and  hence,  it  consists  of  a  number  of  turns  of 
carefully  insulated  wire  wound  into  a  coil  of  known  dimen¬ 
sions.  The  EMP  generated  by  the  sweeping  of  the  magnetic 
field  past  the  probe  is  monitored  on  an  oscilloscope  and  is 
directly  related  to  the  magnetic  field  profile. 

Since  in  cylindrical  geometry  the  self-induceu  field 
is  confined  to  the  azim*'thal  or  9  direction,  the  probe  must 
be  inserted  so  that  its  surface  lies  in  the  p-z  plane  ( see 
Figs.  19  and  27) .  Figures  I9  and  20  are,  respectively,  a 
schematic  of  the  probe  and  probe  circuit;  and  a  photograph  of 
an  actual  B9  probe  ( so-called  from  here  on  because  it  de¬ 
tects  the  magnetic  field  in  the  9  direction) . 

The  pulse  transformer  is  used  so  theit  both  ends  of  the 
probe  may  float.  If  one  end  were  grounded  lureakdown  might 
occur  from  the  inner  cylinder  to  the  probe.  The  relatively 
large  size  of  the  probe  and  the  square  shape  are  the  result 
of  the  probe  theory  presented  in  Appendix  C. 
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3.2.4  Electric  Field  Probes 


As  discussed  in  Chapter  2,  the  radial  electric  field 
ahead  of  the  shock  is  an  important  parameter  in  the  theory 
of  ionizing  waves.  Attempts  were  made  to  measure  this  elec¬ 
tric  field  both  ahead  of  and  behind  the  shock  front.  Limited 
success  was  obtained  in  the  latter  case  using  a  probe  very 
similar  to  the  B0  probe.  The  idea  is  to  insert  two  bare 
wires  into  the  conducting  flow  and  measure  the  potential  dif¬ 
ference  between  them.  Ideally  this  is  possible  because  in 
the  shock  front  and  post-shock  gas  the  conductivity  is  very 
high  and  therefore,  in  the  presence  of  an  electric  field  the 
probe  will  draw  a  small  amount  of  current  and  act,  essentially, 
as  a  voltmeter. 

The  probe  construction  is  shown  in  Fig.  21.  The  exter¬ 
nal  circuitry  is  identical  to  the  B0  probe  (i.e.,  the  Can¬ 
non  connector  goes  to  a  pulse  transformer  and  than  into 
a  length  of  50  ohm  cable  terminated  at  an  oscilloscope) . 

Measurement  of  the  electric  field  ahead  of  the  shock 
front,  i.e.,  in  the  cold,  pre-shock  gas,  was  attengsted  by 
connecting  a  very  high  impedance  voltage  divider  between  in¬ 
ner  and  outer  conductor  and  trying  to  measure  the  potential 
drop  between  them.  As  simple  as  this  may  sound,  the  problems 
are  very  complex.  First,  t.he  voltage  divider  must  be  a  very 
high  impedance  one  to  preclude  any  loading  of  the  shock  tube. 
Second,  since,  before  breakdown,  the  total  bank  voltage  appears 
between  the  cylinders,  the  -voltage  divider  :-.iust  be  able  to 
withstand  upwards  of  50,000  ^'olts  and  alsc  educe  this  50,000 
volts  to  manageable  proportions,  (The  divider  must  divide  a 
few  hundred  to  one,  at  least.)  Third,  s.rnce  the  ...oxrmum  ex¬ 
pected  alectric  field  after  breakdown  is  only  of  the  order  of  a 


>  - 


fev;  tens  of  volts  per  centimeter,  th^  divider,  from  the  second 
consideration,  reduces  this  to  a  very  low  level.  Fourth, 
noise  considerations  require  that  all  measurements  be  done 
in  a  screen  room  20  to  30  feet  from  the  shock  tube.  Hence, 
the  voltage  divider  must  feed  a  cathode-follower  which  can 
drive  terminated  50  ohm  cable  without  too  much  rise  time 
degradation  or  loss  of  signal  strength. 

However,  even  after  recognizing  the  problems,  the  actual 
measurement  indicates  what  appears  to  be  an  excessively  large 
potential  in  comparison  v;ith  predictions  from  the  theory  and 
from  results  of  bare  wire  probes.  The  source  of  this  disparity 
is,  as  yet,  not  understood. 
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3.3  Data  and  Error  Analysis 


3.3.0  Introduction 


For  a  given  axial  magnetic  field  and  pre-shock  pressure, 
comparison  with  the  theory  presented  in  2,2.1  requires  the  si¬ 
multaneous  measurement,  at  a  specific  distance  from  the  knife- 
edge,  of  the  shock  velocity;  the  instantaneous  total  current 
delivered  to  the  shock  tube;  the  self-induced  magnetic  field, 
in  the  0  direction,  both  behind  the  shock  and  behind  the 
rarefaction  fan  that  foJ.iows  the  shock;  the  length  of  the 
homogeneous  gas  sample  l:)etween  the  shock  and  rarefaction  fan; 
and  the  electric  field  ahead  of  and/or  behind  the  shock. 

The  first  two  of  these  quantities  are  separately  measured. 
The  next  three  are  all  determined  from  measurements  (see 

Sec.  3.2.3) >  and  the  last  is  obtained  from  the  electric  field 
probe  (see  Sec,  3.2.4). 


3,3.1  Shock  Velocity  Measurement 


Photomultiplier  outputs  were  monitored  at  the  20,  25, 

35,  50  and  75  centimeter  ports  along  the  shock  tube  length. 

For  a  specific  initial  pressure,  longitudinal  field,  and 
capacitor  bank  voltage,  the  outputs  were  recorded  on  Poleiroid 
film  and  analyzed  in  the  following  manner.  The  photomulti¬ 
plier  trace  at  each  port  was  magnified  and  projected,  using 
an  overhead  projector  v;ith  pre-set  magnification,  onto  a 
single  sheet  of  graph  papex*.  In  this  way  all  five  ports  could 
be  presented  on  a  single  time  base  for  easy  examination.  The 
time  of  arrival  of  the  first  luminous  front  at  each  port  was 
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measured  and  recorded.  Time  zero  was  determined  by  projecting 
the  current  trace  onto  the  same  display  and  using  the  instant 
when  current  first  began  to  flow  as  the  t  =  0  point.  For 
each  set  of  initial  conditions  the  shock  tube  was  fired  an 
average  of  12  times.  The  arrival  times  at  each  port  were 
tabulated  and  plotted  on  an  x-t  diagram  and  a  curve  was 
fitted  to  the  data.  The  slope  of  the  curve  at  any  point  gives 
the  average  shock  velocity  at  that  point.  Figure  22  is  one 
such  typical  x-t  diagram.  Data  is  best  fitted  by  the  heavy 
straight  line  drawn  through  the  points,  while  the  two  dotted 
lines  indicate  the  "steepest”  and  "least  steep"  curves  that 
can  be  sensibly  dravm  through  the  data.  The  shock  velocity  is 
calaculated  as  6.70  X  10®  jt^O^  cm/sec,  and  the  error  is 
assigned  by  measuring  the  slope  of  the  two  dotted  curves.  This 
is  the  procedure  that  was  followed  in  each  series  of  sshots  and 
results  similar  to  Fig.  22  were  always  obtained,  i.e,  between 
20  and  75  centimeters  the  shock  front  v;as  propagating  steadily. 
(Though  we  are  not  sure  the  luminous  front  is  coincident  with 
the  shock  front,  the  fact  that  the  velocity  of  the  luminous 
front  is  steady  enables  us  to  use  that  velocity  as;  the  shock 
velocity. )  In  all  cases  the  error  was  no  more  than  +  12^ 
and  no  less  than  ^  8^.  For  ease  in  subsequent  calculations, 

+  10^  has  been  accorded  each  velocity  measurement.  Figure  23(a) 
is  a  typical  set  of  phocomultiplier  traces  for  a  given  set  of 
conditions.  Figure  23(b)  shows  a  set  of  traces  projected, 
with  the  current  trace,  onto  a  common  time  base. 

Table  I  gives  the  experimentally  measured  shock  velocity 
for  the  series  of  runs  used  as  a  basis  for  this  experimental 
study.  Note  that  the  series  covers  the  sub-Alfvenic,  trans- 
Alfvenic,  and  super-Alfvenic  regions.  The  ordinary  sonic 
Mach  number,  M,  has  been  included  to  give  some  perspective  to 
the  numbers. 
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3.3.^'  Instantaneous  Current  Measurement 


The  instantaneous  current  monitored  by  the  current  moni¬ 
toring  loop  (see  Sec.  3.2.1)  is  displayed  on  an  oscilloscope 
face  and  recorded  on  Polaroid  film.  As  mentioned  above,  it 
is  projected  onto  the  common  time  base  used  for  the  shock  velo¬ 
city  study.  This  serves  a  dual  purpose.  Not  only  does  it 
immediately  provide  a  time  2ero  point  for  shock  arrival  mea¬ 
surements,  but  further  one  can  easily  determine  the  total 
amount  of  current  flowing  through  the  shock  tube  at  the  time 
the  shock  reaches  any  specific  point.  Since  all  magnetic 
probe  studies  (see  Sec.  3.3.3)  were  conducted  at  the  45  centi¬ 
meter  point,  all  the  current  calculations,  now  to  be  discussed, 
give  the  total  current  through  the  shock  tube  for  that  instant 
when  the  shock  had  traversed  45  centimeters. 

Prom  the  x-t  diagram  for  a  specific  run  (consisting 
of  an  average  of  12  shots)  the  most  probable  arrival  time  of 
the  shock  at  45  centimeters  is  determined.  Using  this  time 
the  amplitude  of  each  individual  current  trace,  at  that  time, 
is  tabulated  and  averaged.  With  the  known  conversion  factor 
discussed  in  Appendix  B,  the  average  instantaneot;s  current 
flowing  through  the  shock  tube  is  determined. 

Table  II  gives  the  instantaneous  current  measured  for 
the  runs  listed  in  Table  I.  The  error  associated  with  each 
measurement  is  a  compromise  -ween  the  square  root  of  the 
sum  of  the  squares  of  the  individual  errors  involved  in  read¬ 
ing  a  value  for  l(t)  from  a  ingle  scope  trace  of  finite 
thickness;  and  a  standard  deviation  analysis,  based  on  the  aver¬ 
age  of  all  the  shots,  which  gives  a  somewhat  more  precise 
value.  The  need  for  compromising  lies  in  recognizing  that 
the  current  measurement  ultimately  depends  on  exactly  how  the 
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probe  is  sitting  at  a  given  time,  how  well  the  capacitor  bank 
voltage  can  be  set  and  reset  from  shot  to  shot,  and  other  con¬ 
ditions  not  amenable  to  conventional  error  and  analysis. 

Figure  9  •‘=hcws  a  typical  current  trace  obtained  from  the 
capacitor  bank.  One  mijht  argue  that  its  shape  leaves  much 
to  be  desired,  since  it  is  not  the  desired  square  current 
pulse.  In  fact,  however,  the  ultimate  criterion  for  any  cur¬ 
rent  pulse  shape  in  an  electromagnetic  shock  tube  i s  whether 
or  not  it  gives  steady  shock  propagation.  Since  this  criter¬ 
ion  is  experimentally  well  established  in  our  device  little 
effort  was  made  to  improve  the  pulse  shape. 


3.3-3  Measurement  of  Switch-on  and 
Drive  Magnetic  Field 


In  this  section  we  turn  our  attention  to  the  problem  of 
measuring  and  mapping  the  self-induced  magnetic  field  produced 
by  the  instantaneous  curr-nt  flowing  tlirough  the  electromag¬ 
netic  shock  tube.  The  field  has  two  distinct  components,  each 
of  which  must  be  measured.  The  first  is  the  *'switch-on”  field 
discussed  in  detail  in  Sec.  2.2.1.  It  is  the  magnetic  field, 
pre<licted  from  an  analysis  of  the  ionizing  shock  jump  condi¬ 
tions,  that  appears,  in  the  0  direction,  directly  behind 
the  shock  and  results  from  some  fraction  of  the  total  instan¬ 
taneous  current  flowing  through  the  shock  itself.  The  second 
is  the  *'drive"  magnetic  field  that  results  because  the  remain¬ 
der  of  the  instantaneous  current  must  flow  through  a  simple, 
centered  rarefaction  fan,  i.e.  an  expansive  wave,  traveling 
behind  the  shock.  The  rarefaction  fan  is  not  an  intimate 
part  of  the  problem  of  solving  the  jump  equations,  but  rather 
it  arises  from  the  initial  and  boundary  value  problem  associ- 
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ated  with  the  coaxial  electromagnetic  shock  tube.  The  rare¬ 
faction  fcin  is  needed  to  slow  the  gas  down  behind  the  shock 
so  that  it.  meets  the  boundary  restriction  of  zero  normal  velo- 
;jity  at  the  end  wall,  (Actually^  from  a  mathematical  viev^oint 
one  can  have  either  zero  normal  velocity  or  a  perfect  vacuum 
fit  the  end  wall,  but  the  former  is  probably  the  more  realis¬ 
tic.  ) 


The  mechanics  of  coupling  the  simple,  centered  rarefac¬ 
tion  fan  to  the  shock  jump  problem  and  to  the  boundary  condi¬ 
tions  imposed  by  the  device  are  best  explained  in  Ref.  29  and 
33.  The  former  treats  the  problem  for  ’’classical”  mhd  shocks 
in  a  coaxial  electromagnetic  shock  tube  and  the  latter  extends 
the  analysis  to  the  normal  ionizing  shock  situation.  Here  it 
will  be  sufficient  to  note  that  across  a  rarefaction  fan  all 


fluid  variables  change  continuously,  as  opposed  to  a  shock 
wave  where  one  assumes  a  finite  jump  in  the  fluid  variables. 

The  1  ading  edge  of  the  rarefaction  fan  propagates  into  the 
gas  bihind  the  shock  at  the  characteristic  velocity  associated 
With  that  gas.  If  one  plots  one  of  the  variables  behind  the 
shock,  e.g.  pressure,  at  a  given  time,  as  a  function  of  dis¬ 
tance,  the  profile  appears  as  in  Fig.  24(a).  The  x-t  dia¬ 
gram  for  such  a  situation  is  shown  in  24(b),  the  relation  be¬ 
tween  24(a)  and  24(b)  being  obvious.  In  the  same  situation 
we  may  ask  for  the  current  and  magnetic  field  distribution 
through  shock  and  rarefaction  fan.  The  qualitative  features 
associated  with  each  of  these  quantities  is  shown  in  Fig.  25, 
(b)  and  (c).  Figure  25(a)  is  an  x-t  diagram  showing  a  shock 
of  finite  thickness  and  relates  25(b  and  c)  to  a  specific  shock 
tube  situation.  Note  that  the  total  current  divides  between 
the  shock  and  the  rarefaction'  fan  ana  that  no  current  flows 
in  region  1  chead  of  the  shock,  or  in  region  2  behind 
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the  shock,  or  in  reqion  4  behind  the  rarefaction  fan.  As  a 
result,  the  magnetic  field  in  the  azimuthal  direction  rises 
from  zero  in  region  1  to  a  constant  value  in  region  2,  the 
switch-on  field,  hereafter  denoted  by  Bq( SW) ,  and  then  rises 
again  through  the  rarefaction  fan,  region  3>  a  final  value, 
the  so-called  drive  field  value,  hereafter  denoted  by  BQ(Drive). 
It  is  this  behavior  of  that  we  v/ill  tr^'  and  verify  by 

analyzing  the  output  of  a  probe  (see  Sec.  3.2.3)  located 

at  the  45  centimeter  position. 

In  Appendix  C  a  detailed  analysis  of  the  method  used  in 
obtaining  Bg  from  the  observed  SMF  is  presented.  We  will 
use  the  results  of  that  analysis  for  the  specific  probes  dis¬ 
cussed  in  Sec.  3*2.3-  However,  a  few  words  of  introduction 
are  in  order.  The  essence  of  the  discussion  in  Appendix  C 
is  that  a  large  area  probe  can  be  effectively  utilized  for 
the  point  measurement  of  magnetic  fields  in  moving  media,  for 
the  case  in  which  the  field  variation  is  known  in  two  spatial 
directions  and  unknown  in  the  third.  For  an  infinitely  large 
probe  (effectively,  large  conujared  to  any  spatial  changes  in 
the  magnetic  field)  the  probe  is  equivalent  to  a  single  length 
of  conductor  along  a  direction  of  known  field  variation.  If 
a  magnetic  field  is  moving  across  the  probe,  the  potential 
difference,  or  the  EMF,  between  the  ends  of  the  conductor  is 
(^  X  ^)-t  where  V  is  the  velocity  of  the  moving  field  with 
respect  to  the  probe,  B  is  the  magnetic  field  and  L  is 

the  length  of  the  conductor.  The  EMF  is,  then,  proportional 

dB 

to  B  directly  and  not  to  —  as  is  often  mistakenly  adduced 

^t 

from  a  cursory  examination  of  Faraday’s  law.  For  a  finite 
sized  probe  the  determination  of  B  frcrni  the  observed  EMF 
involves  the  use  of  a  simple  recursion  formula  developed  in 
Appendix  C,  Eq.  Cl^.  An  example  of  the  use  of  that  equation 
will  now  be  reviewed. 
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Assume  the  geometry  depicted  in  Fig.  Cl,  Appen'dix  C, 

and  assume  that  the  SMF  (t)  observed  between  the  terminals 

a  and  b  of  the  probe  is  the  heavy  curve  drawn  in  Fig.  26. 

For  a  given  shock  velocity  U  and  a  probe  dimension,  m,  in 

_  ® 

the  X  direction,  a  recursion  time,  t  ,  rs  defined;  r  =  — . 

Assume  U  and  m  such  that  t  coincides  with  five  divi- 

r 

sions  on  the  time  scale  in  Fig.  26.  The  recursion  formula 

says  one  must  do  the  following.  For  0<Ct<CT,  is  iden- 

r  o 

tical  to  the  EMF.  For  t  <v  t  ^  2i  ,  and  in  oarticalar  for 
point  6,  Bg=[EKF  (at  point  6)  -r  EMF  (at  point  1)],  where 
the  proportionality  constant  is  determined  from  the  probe  di¬ 
mensions  and  the  shock  velocity.  For  2t  \  t  ^ 
particular  for  point  11,  ^0  =  [2^  {at  point  11)  -r  EMF  (at 

point  6)  +  EMF  (at  point  l)j.  In  this  way,  B^  is  mapped  oat 

o  *■ 

point  by  point  directly  from  the  EMF  observed.  The  circled 

points  plotted  above  the  EMF  curve  represent  the  Bg  derived 

from  the  EMF  for  times  t  )>  t  .  The  comolete  B^.  mao  con- 

z'  r  -  6 

sists  of  the  EMF  trace  from  0  -►  t  and  continues  smoothlv 

r 

through  the  circled  points.  The  EMF  shovm  was  constructed 
to  produce  a  B„  such  as  would  be  produced  by  the  electro- 
magnetic  shock  tube  and  shows  all  the  features  we  have  dis¬ 
cussed.  To  convert  the  derived  as  a  function  of  labor a- 

tory  time,  to  Bg  as  a  function  of  its  own  spatial  coordinate 
one  need  only  mult. .ply  the  times  observed  by  the  shock  velocity 
U.  One  further  point  must  be  mentioned.  From  a  purist’s  point 
of  view,  use  of  the  recursion  formula  through  the  rarefaction 
fan  is  not  justified.  Why  it  is  not  is  explained  fully  in 
Appendix  C.  However,  for  any  of  the  experiments  performed  here 
the  correction  to  the  recursion  formula  is  quite  ssall  so  that 
in  all  actual  cases,  described  below,  as  well  as  the  exa.role 
just  discussed,  the  recursion  formula  was  used  to  calculate 
Bg(  SW)  and  Bq( Drive).  The  advantage  in  che  method  just  out- 
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lined  is  that,  in  principle,  since  it  only  involves  counting, 
the  fri^iiehcy  wandWidth  of  the  measuring  system  is  infinite. 
Further,  of  course,  it  eliminates  any  difficulties  that  de¬ 
velop  from  time  integration  of  the  EMP. 

We  now  direct  ourselves  to  the  actual  ej^erimental  situ¬ 
ation.  The  probe  used  to  detect  the  self-induced  magnetic 
field  consists  of  ten  turns  of  Carefully  wound  Formvar  wire. 

The  entire  coil  is  inserted  in  the  annulus  of  the  shock  tube 
and  the  output  is  coupled  through  a  nominally  pulse  trans¬ 
former  to  a  terminated  50  ohm  cable  ( see  Fig.  I9) .  The  exact 
probe  geometry,  size  and  position  is  shown  in  Pig.  27. 

Equation  GI5  of  Appendix  C  states  that  the  proportion¬ 
ality  factor  between  and  the  EWF  is  just  where  U 

is  the  shock  velocity  and  t  is  the  probe  dimensF&n  in  the 
radial  direction.  The  derivation,  however,  assumes  ho  varia¬ 
tion  of  Bci  in  that  direction,  in  the  coaxial  electromagnetic 

”  1 
shock  tube  the  field  obviously  varies  as  —  and  this  can  be 

o  .  p-v 

simply  introduced  into  Eg.  CI5  by  replacing  t  by  (^,n 

( see  Pig  25) .  The  proportionality  Constant  then  becomes^ 

-l/u(-tn  — p  -  and  for  the  probe  shown  in  Fig.  27  we  get 
■  Pi  ‘^mean 

-l/u( i 00163) i  Equation  GI5  for  such  a  one  turn  probe  becomes; 


P _ ,  t)t^^Sf££]  =  — - — y.  EMP 

8  medn  U(.00l63)  ^ 


(15} 


where  Bq(45,  t)  ^  the  self-induced  magnetic  field  in 

the  0  direction  at  time,  t,  at  the  axial  position  45  centi¬ 
meters  from  the  knife-edge,  and  at  the  radial  position, 
from  the  -shock  tube  center  line, 

*  The  summation  term  has  been  left  purposely  vague.  Refer 
to  Appendix  C,  Eg,  G15  for  its  complete  form. 
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For  a  ten  turn  probe  and  a  pulse  transformer  whose  trans¬ 
former  ratio  is  exactly  1:5.2,  the  relationship  between  the 
EMF  observed  across  the  terminated  50  ohra  cable ^  called 


'out 


in  Fig.  19,  and  2^(45,  p, 


mean-* 


t)  becomes: 


Pmean^ 


.52 


U( . 00163) 


'out 


or 


B  p  .  t)iiiS5|££]  .  3,141011  ^ 

S'-  -'■*  »^mean^  ■  lor  ^ 


U[E] 

s 


'out 


(16) 


This  is  the  formula  used  in  all  subsequent  calculations. 

For  each  shock  tube  firing  the  output  of  the  probe 

was  monitored  at  the  45  centimeter  point.  The  polcuroid  record 
was  projected  onto  graph  paper  just  as  were  the  current  aM 
photomultiplier  tracers.  Each  trace  was  reduced  using  the 
analysis  described  above.  The  values  of  Bq(SW),  and 
Bq( Drive),  and  the  length  of  the  hocaogeneOus  gas  sample  were 
recorded.  All  shots  at  the  same  set  of  initial  Conditions, 
conprising  a  run,  as  listed  in  Tables  I  and  II,  were  then  tabu¬ 
lated  and  averaged  to  give  an  average  B^(sw),  B^CDrive), 
and  homogeneous  gas  sample  length.  The  results  are  listed 
in  Table  III,  The  table  is  seif-e:q>lanatory  except  for  the 
last  column,  labeled  .  At  any  time,  t,  after  t  =  Oi 

X  represents  the  distance  from  the  3cnife-edge  of  the  leading 
edge  c£  the  rarefaction  fan,  while  X  represents  the  dis- 
tance  from  the  knife-edge  to  the  shock  front  (see  Fig.  25(a)  ). 
As  is  obvious  from  Fig.  25  (a)  both  X  ai^  X  depend  linearly 
on  t:  and  therefore  the  ratio  X^'^Xg  is  time  ii^epeiKient 
and  represents  a  convenient  way  of  speci^ing  the  length  of 
the  homogeneous  gas  sample  (hgs).  At  ai^  axial  position  along 


(17) 


X  -  X  =  X  (1 - )  . 

c  ^  .Q  XT  * 


The  error  associated  with  each  of  the  measurements  tabu¬ 
lated  in  Table  III  will  now  be  discussed  in  details  Figures 
28,  29  and  30  shov;  three  different,,  but  representative, 
traces.  The  top  trace  in  each  figure  is  the  actual  SMF 
obtained  for  the  conditions  noted,  while  the  bottom  trace 


EMP  and  the  B, 


8 


represents  an  enlargement  of  the  actual 
derived  from  that  EMP(cf  with  Fig.  26).  The  three  derived 
Bq’s  show  distinct  similarities  and  distinct  differences. 

For  convenience,  we  will  refer  to  them  as  Types  28,  29  and 
30,  respectively.  All  types  shov/  a  characteristic  oscilla¬ 
tory  structure,  albeit  in  28,  tha  t  structiire  is  more  pro¬ 
nounced  than  in  t  h  e  o  thers.  Types  28  and  29  both 
show  a  distinct  2'^level  structure,  exactly  as  we  expect  for 
a  switch-on  shock  followed  by  a  rarefaction  fan  (Cf  Fig.  26). 


Type  29  has  an  unambiguous  swdtch-on  regime  while  the  switch- 
on  region  in  28  is  apparent,  but  not  as  dramatic.  In  30 
the  switch-on  region  is  not  well  formed.  We  believe  the  slope 
change  in  the  trace  is  diffusion  of  the  magnetic  field  through 
the  region  between  the  Shock  and  rarefaction  fan.  The  dotted 
lines  in  each  figure  represent  the  average  value  of  the  switch- 
on  and  drive  magnetic  field  for  that  shot.  At  the  present 
time,  it  is  believed  that  the  oscillatory  structure  of  the 
traces  not  part  of  the  physics,  but  rather  the  result  of 
the  electrical  characteristics  of  the  probe  circuit.  In  fact, 
an  investigation  of  the  circuit  points  to  the  pulse  trans¬ 
former  as  the  element  in  the  system  responsible  for  the  oscil¬ 
lations.  To  see  why  this  is  believed  to  be  the  case  consider 
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Fig,  3i-  ^‘ne  heavy  line  represents  an  actual  El^  obtained 
in  a  typical  shot.  The  dotted  curve  represents  the  magnetic 
field  obtained  from  the  EMP,  We  note  the  oscillatory  struc¬ 
ture.  If,  now,  we  draw  a  smooth,  monotoniC  Curve  (the  dashed 
curve)  through  the  derived  Bq  and  re-derive  the  EMF  that 
would  give  this  dashed  curve,  we  get  the  dashed  EMP  curve; 

It  is  apparent  that  some  of  the  structure  necessary  to  smooth 
out  the  oscillations  in  the  derived  from  the  experimental 

EMP  has  been  suppressed.  The  pulse  transformer’s  response^ 
in  fact,  does  fall  off  significantly  at  the  frequencies  repre¬ 
sented  by  the  suppressed  structure  and  for  this  reason  the 
oscillatory  behavior  of  the  derived  magnetic  field  is  not 
thought  to  be  indicative  of  actual  magnetic  field  structure 
in  the  plasma. 

Por  a  given  series  of  shots  the  error  in  the  measure¬ 
ment  of  an  average  Bg(sw)  or  Bg(Drive),  or  X is  a 
combination  of  two  factors.  The  first  is  the  natural  scatter 
that  would  be  obtained  in  any  experiment.  Thus,  even  if  all 
the  Bg(SW)’s  looKed  like  29  would  expact  to  see  a  varia^ 
tion  in  the  measured  level  from  shot  to  shot.  This  kind  of 
Variation  can  be  accounted  for  by  a  calculation  of  the  stai^ard 
deviation  of  the  average  measurement.  The  second  Source  of 
error  is  the  error  associated  with  the  assignment  of  some 
Value  to  an  individual  Bq(sw)  or  B^( Drive)  when  it  is  of 
the  type  28  or  30.  If  the  oscillatory  structure  is  larger 
there  is  a  large  error  in  each  individual  measurement,  dince 
in  any  series  of  shots  there  were  several  of  each  of  the  three 
different  types  of  trace,  it  was  inpossiblei  a  priori,  to  es^ 
tablish  either  error  as  the  dominant  one,  Hence^  the  following 
rule  was  followed.  All  the  ii^ividual  average  \^lues 

(switch-on  ai^  drive)  were  tabulate  for  a  given  series  ax^ 
each  individual  value  was  accorded  its  own  error  to  account 


. . . I . . . 


for  the  oseillatofy  structure.  A  final  average  value  was  cal¬ 
culated,  A  standard  deviation  analysis  was  performed  for  this 


average  value  and  an  error  obtained.  Further,  the  square  root 
of  the  sum  of  the  squares  of  the  individual  uncertainties  was 
computed  and  compared  with  the  standard  deviation.  The  larger 
of  these  two  values  was  taken  as  the  error  in  the  derived  Bg, 
In  the  final  computation  of  the  error  we  must  also  take  into 
account  the  uncertainty  in  the  shock  velocity  since  it  appears 
in  the  final  reduction  of  the  BMP  to  Bg  ( see  Eq.  l6) , 

The  assignment  of  an  error  to  the  measurement 

was  done  only  on  the  '  .sis  of  a  standard  deviation  analysis. 
Examination  of  Figs.  28  and  30  shows  that  it  is  difficult  to 
assign  a  specific  uncertainty  to  each  individual  measurement 
of  the  hgs.  And,  since  for  all  the  cases  just  discussed  the 
standard  deviation  error  turned  out  to  be  the  largest,  our 
neglect  to  perform  an  error  analysis  on  the  individual  errors 
is  justified.  In  addition,  the  length  of  each  hgs  was  so 
small  that  the  natural  scatter  resulted  in  a  large  percentage 
error  in  the  confuted  value  of  X This  is  apparent  in 
the  tabulated  values  presented  in  Table  III  and  this  fact 


relegates  the  parameter  X 
in  the  experimental  study. 


to  a  somewhat  secondary  role 


3.3«^  Electric  Field  Measurements 

As  already  noted  in  3-2.4,  only  limited  success  was  ob¬ 
tained  in  the  measurement  of  the  electric  field  behind  the 
shock  and  almost  no  success  in  the  measurement  of  the  electric 
field  ahead  of  the  shoCk.  The  major  problems  encountered  were 
ij^uctive  pick-up  due  to  the  large  currents  and  short  times  of 
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the  esgjeriment,  and  to  noise  generated  hy  the  air-gap  switch 
firing  circuitry.  Similar  prohleifis  have  heen  encountered 
t>y  other  investigators. 

Various  methods  were  employed  to  eliminate  the  noise 
on  the  probes,  all  without  significant  success.  A  new  approach 
is  being  investigated  at  the  present  time.  For  the  purposes 
of  this  study,  however,  measurements  of  the  electric  field 
must  be  Considered  as  unavailable. 
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3.4  Comparison  with  Theory 


The  theory  of  ionizxng  shocks,  presented  in  2.2.1  leaves 
one  question  completely  unanswered.  In  its  sin^jlest  form,  the 
question  is,  for  any  given  ionizing  shock  velocity,  subject  to 
the  boundary  conditions  imposed  by  an  electromagnetic  shock 
tube,  is  there  a  unique  set  of  dovmstream  variables  consistent 
^ith  ionizing  Shock  theory?  If  we  refer  to  Fig.  3  the  ques¬ 
tion  can  be  rephrased.  For  the  electromagnetic  shock  tube,  is 
there  a  unique  locus  of  points  over  the  range  0  <  u,<  that 
describes  the  propagation  of  ionizing  shocks? 

^h  affirmative  answer  v?ill  help  resolve  the  puzzling 
question  posed  by  the  multiple  contours  characterizing  ioni^ 
zing  shock  theory,  i.e.  which  E  )>  0  Solution,  if  any,  is 
Appropriate  to  a  given  shock  velocity.  This  question  was 
hinted  at,  indirectly,  in  a  discussion  of  the  stability  and 
evolutionarity  of  the  classical,  slow  ihhd  shock  (see  page  20). 

A  single  locus  of  solutions  provides  a  natural  path  of  behavior 
all  ionizing  shocks,  A:^,  because  of  the  general  nature 
of  ionizing  shock  theory,  the  possibility  of  such  a  locus  has 
broad  inm>lications  in  the  many  physical  situations  in  which 
ionizing  shocks  Occur.  For  exan^le,  pinch  ai^  inverse  pinch 
experiments  are  two  areas  of  current  investigation  \diere  ioni¬ 
zing  shocks  play  an  inportant  role, 

^t  will  now  be  shown  tnat  the  esroerrmentai  results  ore— 
ssnted  in  3*  ^  fit  one  such  locus,  and  in  addition,  the  locus 
has  a  distinct  physical  basis  deeply  rooted  in  the  theory  of 
classical  gas  dynamics. 

The  relation  of  ionizing  shocks  to  classical  mhd  shocks 
is  analogous  to  the  relationship  between  detonation  —  deflagra- 


tion  phenomena  and  ordinary  gas  dynamic  shock  phenomena.  In 
detonation  -  deflagration  theory,  the  internal  energy  of  the 
gas  behind  a  reaction  front  (the  detonation  -  deflagration 
analog  of  a  shock)  is  different  from  that  of  the  gas  ahe^  of 
the  reaction  front.  This  difference  in  chemical  composition 
across  the  front  gives  rise  to  the  well  known  ai^  well  veri¬ 
fied  Chapman-Jouguet  hypothesis* The  hypothesis  states  that 
the  entropy  of  the  downstream  gas  assumes  a  stationary  value, 
in  particular,  a  maximum  (minimum)  for  an  expansive  ( con5>ressive 
reaction;  the  gas  behind  the  reaction  front  moves  with  the 
small  disturbance  speed  in  that  gas;  and  the  velocity  of  the 
reaction  front  with  respect  to  the  upstream  gas  also  assumes 
a  stationary  value  with  a  maximum  (minimum)  for  an  expansive 
( coii5>ressive)  reaction.  Further,  the  theory  does  allow  ex¬ 
pansive  reactions,  sinCe  this  is  precisely  idiat  is  n^ant  by 
a  deflagration.  In  ionizing  shock  theory  we  have  seen  that 
the  energy  equation  has  an  electric  field  contfibutipn  added 
to  the  upstreani  energy  in  direct  analc^  with  detpnation-dia^ 
flagration  theory*  The  eff^t  of  this  addition  is,  of  course, 
the  theory  developed  iy  Taussig* 

In  the  context  of  this  discussion  we  esgsand  outline 
of  that  theory  to  incline  the  fellcwing  features*  tn»  icKils 
of  all  et^  points  of  the  S  y  0  contours  las  reaarlg5?le  pr<^i^ 
erties.  For  any  end  jx>iht  of  any  contour,  Taussig  proves; 

(a)  that  the  entropy  of  the  downstream  gas  assi^ss  a 
stationary  value, 

(b)  that  the  velocity  of  the  downstream  gas  is  predsely 
the  slow  disturbance  spe^  in  that  gas,  and 

(c)  that  the  shock  velocity  also  assi^^s  a  staticpal^ 
value. 

These  stat^ents  lead  to  the  conclusion  that  the  locus  of  all 
contour  ei^  ^ints  is  a  locus  of  Chapman-uouguet  ^Ints,  %d.th 


the  following  exception*  The  class  of  end  points  labled  qg 
(Fig.  3)  correspond  to  maximum  entropy  and  maxiraum  shock  velo¬ 
city  points  whereas  the  Chapman- Jouguet  point  in  an  analogous 
reaction  would  have  minima  for  both  quantities.  The  other 
classes  of  end  points,  the  q.^,  and  q^  classes,  satis^ 

all  the  Chapman- Jouguet  restrictions. 

In  Pig*  3  the  dotted  Curves  represent  the  locus  of  the 
E  ^  0  contour  end  points.  In  region  ABC,  E*  is  the  maximum 
value  of  the  electric  field  discussed  in  2,2,1.  The  class  of 
q.g  end  points  all  lie  to  the  “right”  of  E*  while  the  qa 
end  points  ail  lie  to  its  “left*” 

From  the  foregoing  discussion,  it  wx>uld  be  natural  to 
propose  that  the  most  probable  path  of  ionizing  shock  propa¬ 
gation  is  along  the  locus  of  contour  end  points,  starting  at 
U,.  =  0  and  following  the  loCus  through  point  A  ai^  point 
to  point  C*  From  C  we  do  not  expect  the  path  to  back  track 
to  B  and  follow  the  dashed  curve  from  B  for  this>  again, 

poses  the  difficulty  that  for  1  ^  2  there  will  be  mul¬ 

tiple  solution  points,  Ifote  likely,  there  is  a  transition 
curve  from  point  C  to  the  upper  locusj  while  for  u^  ^  2  the 
path  follows  the  upper  en2-point  locus  e^ractly. 

Two  objection^  hoi^ver>  suggest  an  alternative  path  ^diich 
is  a  sli^t  ^^ification  of  this  er^  ^int  proposal.  The 
first  objection  is  tist,  the  contours  characterized  by 

u^  y  i,  the  q.  class  of  eis  ^ints  are  all  entropy  d«::reasing 
ai^  there&re  inadmissible  as  shock  solutions.  This  is  apparent 
frcms  Fig.  3  ^daich  shows  the  entropy  increasing  solutions,  for 

all  to  the  “right®  of  q.,  beyos^i  the  “bracket®  point.  This 
le^s  us  to  the  suggestion  that  ^r  y  i,  the  l<K:as  paral¬ 
lels  the  ei^  point  locus,  but  stays  below  it,  as^  connects 
the  "bracket**  points  on  each  contour.  These  are  points  of 
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no  entropy  change  across  the  shock  and  therefore  are  entropy 
minima  with  respect  to  all  other  admissible  solutions  oh  a 
contour.  The  second  objection  is  that  for  the  q  and  q 

w 

classes  of  end  points  there  is  no  uniform  region  of  hot  gas 
behind  che  Shock,  i.e.  X  /X  =1*  But,  experimentally 

€  5 

X  /X  has  been  measured  and  therefore  the  conjeccure  is  that, 

3 

just  as  iii  the  >  1  Situation,  the  locus  of  solutions  is 
a  path  just  off  the  end  point  locus  and  paralleling  it^  It 
is  well  to  note  that  the  proposal  of  off-ehd  point  behavior 
is  not  unreasonable.  In  early  detonation-deflagration  studies, 
attempts  to  verify  the  Chajmsan-JouqUet  hypothesis  indicated 
definite  off  end  point  behavior.  This  was  later  attributed  to 
I  the  fact  that  the  sound  speed  calculation  used  in  the  theory 

I  was  for  a  single  species  gas  and  not  for  a  mixture  of  species 

I  as  occurs  in  any  detonation  process.  When  multiple  species 

I  gases  were  considered  the  off  end  point  behavior  was  explained. 

I  A  similar  situation  is  probably  responsible  for  such  b^avior 

in  ionizing  shock  studies. 

I 

i  A  basic  aim  of  this  report  is  the  verification  of  the 

proTOsal  that  there  is  a  unique  locus  of  solution  points 
to  the  normal  ionizing  shock  theory  and  that  the  locus  is 
one  that  parallels  the  contour  end  point  locus.  In  so 
doing  we  will  dei^nstrate  uhambigucus  separation  of  the  shock 
I  current  frcsa  the  drive  current  and  the  ^ustence  of  a  small  > 

I  but  distiguishable^  uniform  hot  plasma  sa^le.  This  ^itch- 

I  On  b^avior  of  the  ionizing  front  will  be  shewn  to  extend 

I  to  slew  and  fast  shocks  as  predicted  by  the  normal  ionizii^ 

I  shock  theory . 

I  Table  IV  is  a  cc^^arison  between  experimental  results 

I  and  the  result  of  restricting  the  normal  ionizing  shock 


theory  to  the  prOtKssal  of  off  end  |K>int  locus  b-.uavior.  The 
"General  Statistics**  Section  is  self-expiasatory.  The  ex- 


. . 


perimental  tabulation  has  already  been  discussed.  The  theo¬ 
retical  tabulation  will  be  discussed  now.  Theoretical  value 
were  obtained  as  follows.  For  each  run,  #1  -  11,  the  output 
of  a  computer  program  (designed  to  give  all  downstream  var¬ 
iables  for  a  given  set  of  upstream  variables  and  including 
a  jump  in  7)  was  plotted  over  the  entire  range  of  the  elec¬ 
tric  field,  i.e.  from  E=0  to  the  value  of  E  at  which 
a  contour  end  point  was  reached.  In  reality,  the  exact 
contour  end  point  cam  only  be  obtained  with  extraordinary 
labor  or  extraordinary  luck  since  the  program  employs  a 
successive  approximation  method.  Hence,  the  electric  field 
was  actually  swept  from  0,  not  to  the  end  point,  but  to 
within  1/2%  of  the  end  point,  sufficient  accuracy  for  our 
purposes.  In  addition,  since  the  experimental  error  in  the 
shock  veloc'-ty  was  +  IC^,  each  run  was  supplemented  by 
programs  covering  ±  around  the  central  value.  In  this 
way  the  error  in  the  theoretical  values  could  be  included 
in  the  analysis.  Note  that  each  theoretical  value  in  the 
table  is  followed  by  two  error  figures.  The  upper  one 

corresponds  to  M  +  10^  and  the  lower  to  M  -  IC^.  Tn 

a  a 

some  cases  both  errors  are  of  the  same  sign.  This  will  be 
explained  later  on. 

For  the  first  seven  runs  (runs  -  11  will  be  examined 

separately),  the  following  procedure  was  used  to  obtain 

theoretical  values  and  theoretical  errors.  At  the  computer 

program  end  point,  as  opposed  to  the  exact  end  point,  X  /X 

©  s 

(theoretical)  and  X  /X  (experimental)  were  compared.  In 
all  cases,  they  were  within  3^,  in  most  less  than  1%,  The 
values  of  Bq{  Bq( Drive),  l(t),  and  X  /X  were  then 

read  from  the  plots,  at  the  computer  end  point,  and  recorded 
with  their  respective  errors.  Since  we  propose  that  opera¬ 
tion  is,  indeed,  off  end  point  we  are  within  rights  in  using 
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as  we  did  the  computer  program  end  point.  Ahd>  in  addition,  | 

as  was  verified  in  several  cases,  between  computer  end  point  | 

and  actual  end  point  none  of  the  variables  vary  significantly*  | 

In  several  cases,  one  example  being  in  run  both  | 

errors  are  negative.  This  simply  means  that,  at  this  value  f 

of  Bq(sw}  |theor^  is  a  relative  maximum.  Varying  the  | 

shock  velocity  ±  10^  around  the  value  gives  smaller  switch-  I 

n  fields.  The  reasoning  is  exactly  the  same  for  two  positive  | 

errors,  only  here  the  main  value  is  a  relative  minimum  (see  | 

run  7^,  Bq(sw)  rheory).  When  the  errors  are  of  opposite  f 

sign  but  of  different  magnitude  it  means  that  the  main  re-  § 

suit  is  close  to  a  stationary  value.  1 

In  the  first  seven  runs,  agreement,  though  crude>  is  J 

reasonably  good.  In  only  one  run  (run  is  there  more  i 

than  one  value  that  is  beyond  the  range  of  experimental  arid  M 

theoretical  error.  In  addition,  i**!  runs  #1,  2^  4,  and  5  | 

a  variation  of  less  than  3^  the  off  end  point  position  1 

brings  all  values  within  experimental  error .  .4  disparity  1 

of  3^  is  not  unreason^le  when  viewed  in  the  context  of  ?the  § 

early  detonation-deflagration  studies  Where  a  5^  error  was  ■ 

common.  In  run  a  12^  shift  in  end  point  is  needed  to  M 

bring  Bq(  SW)  into  agreement  and  a  31^  shift  is  needed'  for  j 

B0(sw),  run  #3.  These  last  two  points  are  distictly  put  bf  B 

line.  For  runs  #6  and  Bq (Drive),  in  each  case;,  is  put-  ft 

side  the  bounds  of  experimental  error .  Variatioii  in  the^  | 

end  point  position  cannot  account  for  the  error  and  these  1 

two  individual  measurements  must  be  vi^ed  with  suspicion  *  j 

Runs  #8  and  are  borderline  runs.  In  each  case  i 

M  ,  and  -  10^  are  trans-Alfvenic,  but  M  +  iC^  is  sUppr-  § 

Alfveni-C,  and  in  the  super-Alfvenic  region  the  magnitudes  p 

of  the  downstream  variables  are  quite  different  from  the  p 

S' 

magnitudes  in  the  trans-Alfvenic  region.  This  is  immediately  p 

W 

W: 

SI 


apparent  from  the  errors  assigned  to  the  theoretical  values 
for  the  runs.  The  enormous  range  available  for  the  switch*- 
on  fields  makes  the  measurement  almost  pointless  unless  the 
shock  velocity  error  can  be  reduced  to  less  than  5^* 
any  event  it  is  comforting  that  the  measured  values  of 

both  cases,  lie  above  the  main  theoretical  value, 
as  they  must  if  they  are  to  fit  the  theory  in  any  sense.  The 
drive  fields  on  the  other  hand  are  not  within  the  experimental 
error  and  variation  of  the  end  point  position  cannot  bring 
them  into  line.  In  the  total  picture  runs  and  ^  must  be 
viewed  as  suggestive,  at  best.  Only  the  Bg( Drive)  measure¬ 
ment  does  not  agree  with  experiment,  but  even  so  until  the 
shock  velocity  error  can  be  reduced  the  range  of  switch-on 
error  and  X  /X  error  will  remain  disconcertingly  large. 

Runs  f^lO  and  jSll  are  both  super-Alfvenic  runs.  We  can¬ 
not  go  immediately  to  the  contour  end  points  for  our  theore¬ 
tical  values  because,  as  we  have  seen,  they  are  entropy  de¬ 
creasing  points.  We  must  first  find  the  "bracket”  points 
where  the  entropy  begins  to  go  positive.  However,  in  the 
co.mputer  progfcun  developed  for  the  determination  of  theore¬ 
tical  values,  we  have  used  a  jump  in  7  across  the  shock 
front  and  it  is  not  apparent  how  to  calculate  the  entropy 
jump  for  such  a  situation .  All  that  can  be  proved  from 
examination  of  the  steady  state  shock  theory  is  that,  even 
if  7  jumps  across  the  shock,  the  g^  class  of  end  points 
is  always  entropy  decreasing.  The  exact  numerical  value 
cannot  be  simply  computed.  The  best  that  can  be  done  is  to 
return  to  the  steady  state  theory  where  7  does  not  jump 
and  find  the  value  of  E  at  which  the  bracket  point  occurs . 
With  this  as  a  guide,  a  guess  can  be  made  as  to  where  the 
bracket  point  occurs  for  the  case  where  7  jumps.  In  fact, 
this  procedure  is  perfectly  acceptable  because  the  numerical 
values  of  the  downstream  variables  do  not  differ  by  more 


than  .1^  between  the  case  where  y  is  constant  across  the 
shock  and  the  case  where  7  is  assumed  to  jump  across  the 
shock . 

For  the  runs  under  consideration  it  turns  out  that  the 
bracket  point  in  each  case  is  no  more  than  .5^  from  the 
contour  end  point.  Using  the  value  of  E  appropriate  to 
the  bracket  point  the  theoretical  values  for  Bg(  SW), 

Bp (Drive),  X  /X  ,  and  l(t)  are  listed  in  their  respective 

o  6  S 

columns . 

Agreement  between  theory  and  experiment  in  run  ^10  is 
fair.  Bq( Drive)  and  l(t)  are  within  the  error  bounds  while 
X^/Xg  and  Bg(  sw)  are  not.  However,  a  shift  of  6^  in  the 
bracket  point  will  bring  all  quantities  into  agreement.  In 
run  #1.1  Bq^  )  and  l(t)  are  beyond  the  bounds  of  the 
assigned  errors,  while  Bq( Drive)  auid  are  well 

within  the  tabulated  uncertainties.  A  shift  of  in  the 
bracket  pc.  't  will  make  B0(sw  )  consistent  with  the  other 
results,  but  such  a  shift  merely  accentuates  the  disparity 
in  l(t).  We  are  therefore  forced  to  view  run  #11  as  spine- 
what  anomalous. 
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3 . 5  Conclusions 

Using  a  co-axial  electromagnetic  shock  tube  we  have 
been  able  to  propagate,  record,  and  analyze  steady,  plane 
switch-on  ionizing  shocks  in  the  sub-Alfvenic,  trans-Alfvenic 
and  super-Alfvenic  regimes.  The  results  are  in  confirmation 
of  the  theory  of  normal  ionizing  shocks.  The  experimental 
existence  of  the  super-Alfvenic  switch-on  front  is  strong 
evidence  that  a  recent  argument  suggesting  that  such  shocks 
have  no  structure  may  be  in  errors 

We  have  presented  a  proposal,  modeled  on  the  Chapman- 
Jouquet  theory  of  detonations,  to  describe  the  behavior  of 
such  ionizing  shocks.  It  has  been  shown  that  the  data  ob¬ 
tained  is  in  tentative  agreement  with  such  a  theory,  except 
for  the  region  in  the  immediate  vicinity  of  the  transition 
point  from  trans-Alfvenic  to  super-Alfvenic  propagation. 

Existence  of  a  small,  but  definite,  homogeneous  gas 
sample  has  been  demonstrated  and  shown  to  fit  into  the  pro- 
p  osed  scheme. 

In  conclusion,  the  experiment  described  in  this  report, 
gives  new  insight  into  the  behavior  of  an  important  class 
of  magnetohydrodynamic  phenomena  -  the  normal  ionizing  shock 
front . 

No  attempt  was  made  .o  propagate  expansive  (sub-sonic) 
ionizxng  shocks  even  though  their  existence  is  predicted  by 
the  theory-  It  is  well  known,  from  conventional  shock 
studies,  that  such  shocks  do  not  produce  significant  ioniza¬ 
tion  of  the  post-shock  gas  and  thus,  they  vrould  not  be 
covered  by  the  ionizing  shock  theory.  theory  cannot  be 

expected  to  explain  phenomena  in  this  vexocity  reinge. 
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FIG.  4  THE  MAGNETIC  PISTON  IN  CO-AXIAL  GtOMETRY  . 
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FIG.  8  SHAPED  -  PULSE  CAPACITOR  BANK 


.227x10-  ompt/cm. 


2 .  Representative  Capacitor 
3-  Series  Inductor 


4 .  Enclosure 


FIG.  10  CAPACITOR  BANK  AND  ENCLOSURE 
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TEFLON  BREAKDOWN  GUARD  RING 
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FIG.  1  I  AIR-GAP  SWITCH  AND  FIRING  CIRCUIT. 


3-  Teflon  Breakdown  Guard  Ring 
Teaser 
5-  Cat5K>de 
6.  Rlexiglas  Support. 

7-  Main  Body  of  Sfe>ck  Tube 
8.  Current  iksnltor  Probe  (See  Fig.  15) 
9-  Ctold  CatiMsde  Gauge 
10,  Magnet  C^ils 


FI6.  12  AIR  GAP  SWITCH  IN  OPERATING  POSITION  1 
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F!G.  13  30  Me /s  ,  R.F.  PRE- IONIZER  CIRCUIT. 


86 


FIG.  14  TIMING  AND  FIRING  CIRCUITS. 
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FIG.  18  PHOTOMULTIPLIER  CIRCUIT 
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FIG.  19  B6  probe  :  schematic  AND  CIRCUIT. 
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FIG,  20  B0  PROBE:  CONSTRUCTION 
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nC.  21  ELECTRIC  FIELD  PROSE:  CONSTRUCTION 


SHOCK  FRONT  ARRIVAL  TIME  IN  fiSEC. 
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FIG.  22  X-t  DIAGRAM  FOR  DETERMINATION  OF 

SHOCK  VELOCITY  . 
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(b)  CURRENT  AND  PHOTOMULTIPLIERS  PROiECTEO  ON  A 

COMMON  TIME  BASE 


FIG.  23  ANALYSIS  OF  PHOTOMULTIPLIER  DATA 
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(b)  CURRENT  DISTRIBUTION  THROUGH  SHOCK  AND 
RAREFACTION  FAN  AT  t| 


(c}  SWITCH  >ON  AND  DRIVE  MAGNETIC  FIELD  AT  TIME  ti 
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FIG.26  SAMPLE  EMF  a  Bq  DERIVED  FROM  IT 
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FIG. 27  SINGLE  TURN  Bg  PROBE  GEOMETRY. 
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FIG.  28  Be  analysis:  Bxi  =4,000 gauss, 
Pi  =  ZOOft,  Ug  =  5.87  x  10®  Cm/Sec 
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FIG.  30  Bq  analysis:  Bxi  =  I.SOOgouss, 
p,  =  200^,  Us=  7.80x  10®  Cm/Sec 


FIG.3I  Bfl  STRUCTURE  STUDY 
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7.  Appendices 


7.i  appendix  a 

ANALYSIS  OF  THE  CAPACITOR-INDUCTOR  BANK 
WHEN  USED  AS  A  SQUARE-PDLSE  CURRENT  GENERATOR 

To  gmerate  an  MHD  shock  wave  in  a  co-axial  shock  tube, 
it  is  desirable  to  have  a  source  of  electrical  energy  capable 
of  producing  a  square,  high-current  pulse  of  10-2C  lisec  dura¬ 
tion. 

The  sin® lest  arrangement  consists  of  capacitors  and  in¬ 
ductors  arranged  i:a  a  low-pass  filter  geometry  (Fig.  Al) . 

This  set-up  simulates  an  open-er^ed  transmission  line 
with  an  input-in® edance  of  fh  and  a  characteristic  discharge 
time  T  =  Zn/LC. 

To  see  hew  this  design  is  related  to  the  desired  square- 
pulse  shape,  we  f  irst  investigate  the  transi€:nt  proper¬ 
ties  of  transmission  lines  and  the  properties  of  the  pulse- 
forming  netvforks  derived  from  them. 

The  usual  approach  to  this  problem  is  from  the  Laplace 
transform  viewpoint. 

Consider  first  the  circuit  in  Fig.  A2.  Kirchoff’s  Law 
requires  that  at  any  time  after  the  switch,  S,  is  closed 


(Al) 
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90 


Multiplying  (l)  by  and  /(  )dt  we  get 


Z  J  i(t)e~^^  +  Rt  /  i(t)e  ^^dt  +  -i-  /  < /  idt  e**°^dt  =  0.  {  A2) 
o  b  %  ^  ^ 


Defining:  i(p)  =  /  i{t)e~°^dt,  we  get 

o 


oo  z. 


i(p)Rj^  +  i(p)Z^(p)  +  ^  /  /  e  ®^dt  =  0  . 


(  A3) 


N  o 


From  Goldman^  Transf.  Calc,  and  Elec.  Trans..  Prentice-Hall, 
Inc.,  19^9 j  pg.  61,  we  have 


7%dt^  e-P^dt  =  ^  +  fholl 

o  s  s 


where 


00 


F(s)  =  /  f(t}e"®^dt 


f  ~(t)  =  /  f{t)dt. 


We  now  have  (3)  in  the  form 


««!■ 


+  Z„  +  -i^)i(p) 


N 


=  -  si 


0^1 


C  o' 

jr 


-sP 


( A4: 


where  q{o^)  is  the  charge  on  the  capacitor  before  the 

switch  is  closed.  Since  -  -  -V„  =  initial  ootential 

droD  across  the  capacitor  we  can  write  (4)  as 


where 


'"o  =  -Vs 


A3; 


i 
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and  2^  is  the  network  impedance  (which  v;iil  be  a  function. 
as  yet  unknown,  of  the  Laplace  parameter  p) . 

Let  us  now  require  that  the  discharge  current-pulse  shape 
be  a  single  square-pulse  of  duration  t  and  magnitude  I,  i.e., 

j (t)  =  I  0  ^  t  \  T 

i(t)=0  t^O;  T^t'Coo. 

irom  the  definition  of  the  Laplace  transform  for  i(t) 
we  get 


i(p)  = 


Substituting  this  expression  into  (5)  and  solving  for  ( ) 

V 

-1  +e-E") 

-  (a6) 


z  + 


(1  -e-P^) 


which  can  be  put  into  the  form 


C 


1  _ 


nP 


=  R. 


V  ^  , 

(nr- 

coth  HL  +  ^ 

2 


EL  zEI. 

b  a  _  2— 


EL 


(A?) 


by  multiplying  (6)  through  by  e  2  in  the  num  .ator  and  de¬ 
nominator. 

If  we  now  choose  =  2IRj  then 


^  cfe  =  h,  ^ 

N 


(A8) 


which  is  just  the  Laplace  impedance  function  for  an  open-ended 
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lossless  transmission  line  of  characteristic  impedance  =  «L' 

and  one  way  transmission  time  6  =  — . 

2 

Thus  if  we  want  a  square  voltage  (or  current)  pulse  of  dura¬ 
tion  T  developed  across  the  resistor  we  can  replace 

Zjj  +  — L_  by  an  open-ended  transmission  line  with 

CjjP 

utilizing  this  result  we  return  to  (5)  and  substitute 

”  ^o  hereby  hoping  to  investigate  the  effects 

Kr 

Of  a  transmission  line  discharging  through  a  resistor  where 
may  or  may  not  match 
From  (5)  we  solve  now  for  i(p)  and  get 


i(p)  = 


p(Rj^  +  coth  p5) 


(  A9) 


Expanding  the  expression, 

V  '  /^  -  -  /Z  -  R,'\^  .  . 

i(p)  =  _ 1-  -£ - 2;.e-2p6  +  (_o - il'l 

+  *4,)  L  '-’"p  *  +  V 


...  i  ( AlO) 


and  inverting  this  (from  tables) 


i(t)  =  (l-U(t-26)  -  [ ■— £-]  [ U(  t-26 ) -U( t-46 )  j  + 

o  ^  o  ^ 


[ ^ j ^ [ U( t-4? ) -u( t-66 )  j  -  ... 


(  All) 


where 


U(At)  =  1 


for  At  0 


U(At)  =  0 


for  At  <(  0  ana  At  =  (t-n6) 


For  we  get  contributions  to  i(tj  from  only  the 

first  two  terras  of  (ll)  and 

i(t)  =  ^  i“U(t-25)^>  (a12) 

2Rl  . 

which  is  a  simple  rectangular  pulse  of  duration  t  =  26.  Mis¬ 
matching  the  load  gives  rise  to  contributions  from  succeeding 
terms  in  (li)  and  therefore,  effectively,  we  get  reflections 
at  the  load  ( see  Glascoe  and  Lebacqz,  Pulse  Generators.  McGraw- 
Hill,  1948,  p.  178). 

Recognizing  the  advantages  inherent  in  the  transmission 
line  approach  we  inquire  if  v;e  can  design  a  line  to  generate  a 
20  M-sec  pulse.  For  t  =  20  M-sec,  5  =  L  =  transmssion  time  = 

10  iisec.  For  a  representative  signal  velocity  on  the  line  of 
500  ft/4,sec  we  would  need  a  line  of  length  L  =  5OO  ft/jisec 
X  10  M-sec  =  5000  ft,  an  unweildy  laboratory  dimension. 

It  is  apparent  that  we  can  only  shorten  the  line  length 
needed  if  we  use  a  lumped -parameter  network,  where  we  get  high- 
voltage  stability  and  can  use  large  circuit  parameters.  In 
this  process  we  will  lose  our  precision  square  pulse  charac¬ 
teristic  and  the  departure  from  "squareness"  will  be  larger 
the  more  we  make  our  network  depart  from  the  distributed  para¬ 
meter  model. 

Our  discussion  will  now  be  centered  upon  analyzing  the 
network  of  Fig,  AL  If  we  connect  our  "bank"  to  a  battery  we 
have  Fig,  A3. 

Solving  this  for  the  mesh  currents,  we  get  n  difference 
transform  equations,  the  general  one  (except  first  and  last) 
being  of  the  form: 


The  general  solution  to  this  equation  is 

ij-(p)  =  Ae^‘^  +  Be~'^® 

where  A  and  B  are  arbitrary  constants, 
into  (13)  yields  the  condition 

cosh  0  =  1+  (^)p“  . 


(  Al4) 


Substitution  back 


(  A15) 


If  v>e  solve  for  A  and  B  by  considering  the  first  and  last 
meshj  we  can  derlvs  the  input  impedance  transform  for  the  net¬ 
work  which  is 


Zi  (P,n)  =  _ =  JL  [g-inh  (r.  +  1)0 

pii(p)  sinh  (n0) 


{  Al6) 


If  we  now  include  a  resistor  R_  =  /IT  in  the  circuit  of 

^  V  c 

Fig,  C  between  the  battery  and  the  first  inductor  we  will  have 
the  general  case  of  a  lumped-parameter  line  terminated  in  a 
matching  impedance.  The  input  impedance  transform  for  such 
a  circuit  is  simply  Ik  ^  ^  .(p  n)  immediately 

t 

from  the  form  of  (l6)  that  in  this  case 


ii(p)  = 


Vp 


/L  ^  r  sinh  (n  +  l)0 
C  Cp  sinri  n0 


(  A17) 


Pn/LC  +  [ 


But  from  cosh  0  = 


1  +  {LC)p?  we  have 
2 

p  =  — ^  sinh  — 


\  Al3  ) 
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and  substituting  this  into  ( I7)  W"  get  the  general  current 
transform  through  a  ma.chei  load  conn^^cted  to  a  lumped  para¬ 
meter  transmission  line. 


(p)  =  sinh  n6 


sinh(n  +  1)6  +  (2  sinh  ^  -  l)  sinh  n6 

2 


{ai9) 


When  inverted  this  gives  i.(tj  as  the  desired  result.  Hcv?- 
ever,  in  general,  this  is  difficult  to  do  and  usually  requires 
machine  calculations.  When  n  =  5  and  R  =  li  this  has 

L  c 

been  done  and  is  discussed  by  Glascoe  and  Lebacqz,  op.  cit. 
p.  183«  along  with  oscilloscope  traces  showing  the  square 
character  of  the  charging  pulse. 

Although  n  =  p  does  not  give  a  particularly  good  pulse 
shape  the  principle  when  extended  to  larger  n  will  give  in¬ 
creasingly  better  results.  If  now,  instead  of  restricting  the 
bank  geometry  to  the  low  pass  filter  form,  one  looks  for  sys¬ 
tems  that  give  overall  better  performance  in  terms  of  rise¬ 
time  and  flatness,  one  finds  that  the  Type  C  Guillemin  network 
strikes  an  excellent  balance  between  performance,  simplicity 
and  ease  of  assembly.  Without  belaboring  the  point  (all  the 
details  can  be  found  in  Glascoe  and  Lebacqz),  it  is  sufficient 
to  note  that  the  capacitors  in  the  four  meshes  are  not  identi¬ 
cal.  As  Fig.  (  S  )  shows,  ^  )>  )>  C^.  Physically 

this  means  that  when  a  switch  between  bank  and  load  is  ’’closed" 
most  of  the  energy  in  discharges  immediately,  giving  a 

sharp  rise  to  the  pulse,  then  C3,  C^,  and  C,  are  used  to 
provide  the  flat-top  and  high  current  capacity. 
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FIG.Al  LOW-PASS  FILTER  TYPE  PULSE  NETWORK. 
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7.2 


APPENDIX  B 


THEORY  OF  OPERATION  OF  THE  CURRENT  MONITOR  PROBE 

In  the  cylindrical  geometry  shown  in  Fig.  { l6  )  we  have 
from  Ampere’s  Law  (using  m.k.s.  units). 


Cylindrical  -syinmetry  implies 


or 


27rp  =  I(t) 
^o 


3(p,t)  =  -2 - 

27tp 


=  2  X  10 


-7  I{t) 


(Bl) 


(B2) 


(B3) 


The  EMF(t)  generated  in  the  loop  due  to  the  time  variation 
of  B(p,t)  is 


EMF(t)  =  -  A-  /  B{p,t)‘dA 


at 


(54) 


loop 


From  ( 3 ) 


EMF(t)  =  - 


dl(t}  P2 

(2  X  10-^)  J  ^  dp  /  dz 


dt 


Pi 


z. 


Whence 


F.MF(t)  =  -  Piili.  (2  X  10“'^)tn  —  ( 


-  z,  ) 


at 


Pi 
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From  Fig.  (i6 ) 


p2  Pp 

—  =  1.25C;  In  —  =  .22 

Pi  Pi 

z  -  z,  =  10"^  meters. 

^  1 


Therefore 

E?4F(t)  =  -  ^  (2  X  10-')(  .22)(i0-^)  .  (b6} 

dt 

At  the  inpur  to  the  operational  integrator  (6)  is  reduced  by 
a  factor  of  10  due  to  the  10:1  attenuator. 

Therefore 


e.  (t)  =  - 
in'  ^ 


^  (.44)10“^° 
dt 


The  output  of  the  integrator  is 


^out 


(•-.) 


i  /  e^„(t)dt 
1  f 


{B7) 


(B8) 


where  is  the  input  resistance  to  the  integrator  and 

is  the  feedback  capacitance  around  the  operational  amplifier. 
For  =  lO'^^fi  and  =  .OOlixf  (values  appropriate  for  in¬ 
tegration  of  approximately  10  as  pulses)  we  get 


"out 


(t)  =  - 


^  .44)10"-°dt 
dt 


=  .44  X  10“^  I(t) 


(B9) 


Therefore:  l(t)  -  -  ■  x  10”®  e  .  =  .227  x  10®  e  , 

'  '  j|.  ij.  out  ‘  out 


or  finally, 

l{  t) ( amperes) 


2.27  X  io={2I£|£®£)e  (volts) 
volt  out. 
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7.3  appendix  c 


MEASUREMENT  OF  MAGNETIC  FIELDS  IN  MOVING  MEDIA 


A  fundamental  physical  quantity  in  many  plasma  studies 
is  the  self-induced  magnetic  field  produced  by  plasma  currents. 
The  proper  interpretation  of  the  voltage  generated  xn  a  simple 
coil  probe  by  a  magnetic  field  contained  in  a  plasma  wliich  is 
moving  with  respect  to  the  probe  is  often  complex.  An  impor¬ 
tant  physical  example,  the  use  of  magnetic  field  probes  in  an 
electromagnetically  driven  shock  tube,  is  examined  here. 

In  a  shock  produced  plasma  with  c  0  the  magnetic 

field  is  assumed  to  be  frozen  in  the  moving  hot  gas  and  the 
problem  of  mapping  the  magnetic  field’s  structure  is  reduced 
to  the  problem  of  correctly  interpreting  the  effects  produced 
by  a  magnetic  field  sweeping  past  a  stationary  monitoring  in¬ 
strument.  The  central  idea  behind  the  measurement  is  Faraday’s 
Law,  i.e.  to  measure  the  EM?  induced  in  a  coil  of  known  dimen¬ 
sions  due  to  the  changing  magnetic  field.  Once  obtained,  this 
EMF  must  be  related  to  the  magnetic  field  structure  in  the 
plasma.  However,  the  re’.ation  of  the  observed  EMF  (v;hich  is 
a  function  of  time)  to  the  actual  magnetic  field  requires 
careful  analysis  of  Faraday’s  Law  in  moving  media. 

There  is  a  general  tendency  to  oversimplify  the  problem 
by  incorrectly  reasoning  that  since  Faraday’s  Law  predicts  an 
EMF  proportional  to  —  straightforward  integration  of  the 
EMF  will  yield  !j>{t)  °^and  hence  B(t),  since  the  area  of  the 
pick-up  coil  is  fixed.  The  following  discussion  should  help 
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temper  this  tendency  by  showing  specifically  under  what  con¬ 
ditions  such  time  integration  is  valid. 

Consider  a  magnetic  field  frozen  into  a  high  conductivity 
plasma  which  is  moving  with  a  constant  velocity  Ui  in  the 
-f  X  direction  with  respect  to  the  laboratory  coordinate  sys¬ 
tem.  A  prime  on  any  quantity  means  that  quantity  is  to  be  con¬ 
sidered  in  the  frame  fixed  in  the  moving  plasma  while  unprimed 
quantities  refer  to  the  lab  frame.  Let  the  magnetic  field  be 
directed  in  the  -5-  z’  direction  (or.  what  is  equivalent,  the 
T  z  direction  since  the  relative  motion  between  the  frames  is 
only  in  the  x  direction)  and  let  it  be  a  function  of  x- 
only. 

Consider  a  one  turn  rectangular  coil  of  x  dimension  m 
and  y  dimension  (i.e.  Area  =  Im)  at  z  =  0  with  one  leg 
at  X  =  0  (call  this  the  left  leg)  and  the  other  at  x  =  m 
(call  this  the  right  leg).  This  situation  is  depicted  in  Fig. 
Cl. 

As  the  magnetic  field  B{x’)  moves  across  the  coil  there 
will  be  an  EMF  generated  between  the  probe  coil  ends,  a  ard 
b.  Kow  can  we  relate  this  observed  EMF  to  B(x*)? 

In  the  unprimed  frame  we  have  from  Faraday's  Law 

SMF(t')  =  -  —  //  1  ^  (Cl) 

dt 

which,  applied  to  the  coil  geometry  depicted  in  Fig.  Cl  implies, 

d  ^ 

SMF(t)  =  -  —  /  J  B(x,t)  d  X  dy  {C2) 

x=0  y=G 

The  integrand  in  1C2)  Isa  function  of  both  x  and  t  (unprimed) 
because  of  the  relative  motion  of  the  two  reference  frames. 
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Since  B(x, t)  is  independent  of  y  we  get. 


/=} 

EMF(t)  =  -  I  —  f  B(x,t)  dx  (C3) 

dt  k=0 


and  since  the  limits  of  integration  in  (C3)  are  time  indepen¬ 
dent,  we  can  write 


=  G.  Using 


But  by  the  original  frozen  flux  assumption* 

this  and  noting  that  - =  —  we  get. 

3x>  3x 


W  ^  4 


3t  ’ 


EMF(t)  =  -U  t 


m 

/ 

x=0 


dx 


Integration  gives 


EMF(t) 


~U  t  [B(0,t)  -B(ra,t)] 


(C9) 


(CIO) 


Equation  (CIO)  shows,  immediately,  that  even  for  this  idealized 
case  the  EMF(t)  is  not,  in  general,  proportional  to  the  time 
derivative  of  B,  but  is  rather  proportional  to  the  difference 
between  the  field  at  the  left  and  the  field  at  the  right  hand 
side  of  the  coil.  It  is  only  in  the  appropriate  small  probe 
limit  that  one  can  validly  approximate  (C9)  by  a  time  deriva¬ 
tive  of  B.  An  exact  discussion  of  this  is  deferred  until 
later. 

We  now  inquire  as  to  how  Eq.  (lO)  can  be  used  to  map 
out  the  magnetic  field  in  the  plasma.  The  key  to  this  lies 
in  noting  that^  since  B  is  a  function  of  x’  only,  the  field 
maintains  its  profile  in  moving  across  the  coil. 

Thus,  the  field  at  time  t^  that  exists  at  x  ---  0  will 
be  at  X  -  m  at 


*  The  physical  conditions  required  for  frozen  flux  in  shock 
produced  plasmas  have  been  examined  by  A.  Kantrowitz,  AVCO 
AERL  RR  l4l  Oct.  1962.  Also,  Ipth  Solvay  Conf,  Belgium  1962, 
Essentially,  the  magnetic  Reynolds  number  should  be  large 
con^ared  to  one. 


For  simplicity  assume  that  t  =  0  when 
impinges  on  the  left  leg  of  the  probe.  Then, 


just 


0  <  t,  <  H 

u 


EMP(t^)  =  -U  tB(0,t^) 


(C12) 


since  B(m, t,  )  =  0  until  t,  =  —  when  B(x')  first  reaches 

U 

X  =  m. 

For  _  ^  t2  <  —  , 

U  U 


EMF(ta)  =  -U  -B(m,t^)] 


(C13) 


But  from  (cll) 


B(m,t2)  =  BfO.t^  -  5  =  B(0,ti)  . 


Therefore, 


EMFCt^)  =  -U  l.[B(0,t2)  -B(0,t^)] 


(C14) 


But,  from  (C12) 


-U  t  B(  0,  tj^ )  =  EMF(  tj^ )  whence  we 


EMF(tJ  +  EMF(t,)  =  -U  t  B(0,t^). 


Generalizing, 


^EMF  (tj^) 


B(0,tj  = 


(CIS) 


where  t.  is  related  to  t,  by 

i  - 


t.  =  t,  +(i  -  1)(S) 
i  -  U 


i  =  1,  . . .  n 


(C16) 
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and 


0  <  t,  <  2 


Equation  (CI5)  allows  us  to  map  out  b(0, t)  for  all  times 
t  )>  0.  In  practice,  the  EMF(t)  and  U  are  observed  experi¬ 
mentally,  while  m  and  t  are  known  dimensions  of  the  probe. 


Before  discussing  the  refinements  of  the  above  theoiy 
necessary  for  the  interpretation  of  magnetic  probes  in  situ¬ 
ations  where  B  =  B(x’,t’),  let  us  return  to  the  question  of 
whether  one  is  ever  justified  in  simply  integrating  the  EMF(t) 
to  get  b(0, t).  Assume  that  B(x*)  rises  from  an  initial 
value  of  0  to  some  constant  value  B^  over  some  distance 

An 

measured  in  its  own  frame.  If  ~  is  essentially  constant 

ox  ,  . 

over  the  probe  width  we  may  approximate  Eq.  (C9)  by. 


EMP(t)  =  +  Ut(~}  (  /  dx 


Sx' 


r 


x=(r  0 


=  +  Utm{~) 
dx 


x=0 


=  +  U(^)  A 
ox 

x=0 


(cir) 


where  A  is  the  area  of  the  coil. 


But  we  have  already  seen  that 

-U  ~  — 

6t 


whence 


EMF(t)  =  -A  —  (0,t) 
ot 


or 


B(0,t)  =  -  £  /  EHF(t)  dt 
A 


(CI8) 

(C19) 
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The  error  involved  in  the  approximation  (CI7)  is 

.a-B 

ax^f 


.a^B 

ax2 


where 


max 


is  the  maximum  value  of  the  absolute  value 


max 


of  the  second  spatial  derivative  of  B,  and 


m 


is,  as  before. 


the  probe  width.  Thus  the  small  probe  limit  is  only  justified 
when  <«  1  where  a  is  the  distance  over  which  ^  has 

®  ax^ 

a  significant  value. 


The  above  discussion  has  centered  on  the  situation  where 
one  can  transform  from  the  lab  frame  to  a  reference  frame  in 
which  the  magnetic  fie3d  '.s  a  function  of  x'  only.  In  a 
magnetically-driven  sho-:h  tube,  however,  there  is  a  region  in 
which  the  magnetic  field  is  distributed  throughout  a  rarefaction 
fan  and  consequently  there  is  no  reference  frame  in  which  B 
is  totally  time  independent.  An  analysis  will  illustrate  how 
one  can  recover  B(x, t),  when  we  have  a  time  dependence  in 
the  fluid  frame.  This  is  done  for  the  special,  but  important, 
case  of  a  rarefaction  fan  which  is  a  simple  centered  wave. 


The  x-t  diagram  for  an  electroraagnetically-driven 
shock  tube  is  shown  in  Fig.  C2,  where  region  O  is  the  cold 
upstream  gas;  region  1  the  shock  heated  gas  where  B  is  a  con¬ 
stant  =  B, ;  region  2  is  the  simple  centered  rarefaction  wave 
where  B  =  B(— )  ;  region  3  is  the  constant  state  behind  the 

wave  where  B  is  aaain  constant  =  B  ;  where  x,  and 

max 

Xj^  +  m  are  the  left  and  right  hand  sides  of  the  pick-up  probe, 
re^ectively;  and  where  U,  U?  and  are,  respectively,  the 
shock  velocity,  the  velocity  of  the  leading  edge  of  the  rare¬ 
faction  wave,  and  the  velocity  of  the  trailing  edge  of  the 
rarefaction  wave. 

to  time  t^,  (see  Fig.  C3)  the  probe  sees  a  magnetic 
field  of  the  time  independent  variety  discussed  above,  and 
hence  the  above  analysis  applies  applies  up  to  ti ,  After 
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however,  the  probe  experiences  the  time-varying  field  in 
the  rarefaction  wave  and  we  must,  therefore,  re-examine  the 
effects  produced. 

Referring  to  Fig.  C3,  which  is  a  detail  of  Fig.  C2,  we 
see  that  for 

tj.  <  t  <  t2 

the  left  hand  side  of  the  probe  is  in  the  rarefaction  fan  while 
the  right  hand  side  is  still  in  the  time -independent  region. 

Thus,  the  rrght-hand  side  generates  a  voltage  =  while 

the  left-hand  side  generates  a  voltage  =  B{x^,t)t  where 
<  Ux,  U2  is  the  velocity  of  the  characteristic  catting 
Xi  at  time  t,  <  t  <  t2.  The  EHF  actually  observed  is  the 
difference  between  these  two  voltages  since  each  generates 
the  same  polarity  voltage  and,  therefore,  we  have 

EMF(ti  <  t  <  ta)  =  Ux,  B(x-,t)  t  -UB-^(m,t)t  (C20). 

Assuming  we  know  U  ,  we  know  B(x,  ,t,  <  t  because 

we  know  B, ,  t,  and  U.  That  U„  is,  in  fact,  known  remains 
to  be  discussed. 

Continuing  to  the  region  t2  <  t  t,  and  realizing 
that  B  =  that  B(x,  ,  t,  <  t  <  t^)  is  tran^or- 

ted  to  (Xj^  +  m)  and  appears  there  for  "C  t  \  t^.  Knovnng 
B{x^  t  ra,  t2  <  t  <  t3)  we  calculate  b{xt,  t2  ^  t  <  t2)  as 
follows; 

EMF{t2  <  t  <  t^)  =  Uxjt)  B{X;,t)t  -  3{x3^-rm,t)^  {C2l) 

where  U  (i  =  1,2)  means  the  characteristic  velocity  in 
^i 

the  rarefaction  wave  at  time  t  at  position  x^,  (Note; 

X2  =  -r  m) . 
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But, 

+  m,  <  t  <  t3)  =  B(x^,  <  t‘  <  tg) 

where  B(x,  ,  t'  t^)  represents  the  already  calculated 

value  of  B  at  and  at  time  t^^  ^  t'  K  ^2  that  is  con- 

vected  to  +  m  at  time  t^  <C  t  <  tg.  Consequently,  we 
get 

EMF(t^  <  t  <  t3)  =  (t)  B(x,,t)t  -U^  (t)  B(Xi,tx  <  t‘  <  t.) 

(C22). 

In  this  way  B(x,  ,t)  is  calculated  for  all  t  )►  t,  . 

There  remains  for  us  to  examine  how  one  determines  . 

Taussig^^  has  shown  that  for  any  reasonable  situation  in  an 
electroraagnetically-driven  shock  tube  U  decreases  almost 

linearly  between  and  Ug.  If  we  know  and  U,  we 

linearly  interpolate  at  any  t  between  U,  and  and  get 

Uy  at  that  time.  To  determine  U,  and  U  we  note  that 

-  2 

for  the  situation  depicted,  the  field  starts  at  O  and  even¬ 
tually  for  t  t,  reaches  a  constant  value  3  .  Therefore, 

beyond  t.  there  is  no  SMF  generated  across  the  probe. 

In  the  lab  we  see  the  EMP  rise  to  a  maximum  and  then 
fall  to  0  at  t^.  If  one  measures  the  velocity  of  this  point 
at  which  the  EMF  goes  to  0  one  has  measured  U^.  To  get  U, 
map  B{x,,t)  using  (CI5)  for  t  >  0  until  the  point  is 
reached  where  B,  begins  to  depart  from  irs  constant  value. 

This  is  the  point  that  corresponds  to  the  leading  edge  of  the 
rarefaction  fan  and  the  velocity  of  this  point  is  U-  .  Thus 
one  gets  U,  and  for  use  in  measuring  B  for  all  t. 
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FIG.  Cl  PROSE  GEOMETRY  AND  MOVING 
FIELD  CONFIGURATION 


-1^ 


FIG.C3  DETAILS  OF  x-t  DIAGRAM 


7.4  APPENDIX  D 


SHOCK  WIDTH  i\ND  SHOCK  FORMATION  STUDY 

As  a  by-product  of  the  study  described  in  the  text,  we 
present  here  two  interesting  results. 

The  first  concerns  an  estimate  of  the  width  of  the  ioni¬ 
zing  shock  front  and  rarefaction  fan.  From  analysis  of  the 
Bq  data,  we  conclude  that  in  all  cases  the  shock  width  is  of 
the  order  of  2  cm.,  while  the  rarefaction  fan  is  of  the  order 
of  4  -  5  cni.  From  the  data  of  Table  IV,  we  can  derive  that 
the  hgs  is  of  the  order  of  2  -  3 

The  second  result  involves  a  study  of  ionizing  shock  for¬ 
mation  time.  From  gas  dynamics,  it  has  been  shown®^  that  the 
stronger  the  shock  is,  i.e.  the  higher  the  ratio  of  shock  speed 
to  sound  speed,  the  more  rapidly  it  forms;  i.e,  it  takes  fewer 
mean  collision  times  to  reach  a  steady  state.  We  would  expect 
to  observe  this  same  general  feature  in  ionizing  shocks,  and 
this  is,  indeed,  the  case.  Fig,  D1  shows  a  composite  study 
of  three  different  shocks,  ranging  from  weak  (M  ]>  .2)  to 
strong  (M  =2.3).  It  should  be  noted  that  here  the  relevant 

cl 

ratio  is  wave  speed  to  Alfven  speed.  The  strong  shock  is  well- 
formed  by  20  cm.,  while  the  weak  shock  is  barely  formed  at  35  cm. 
The  intermediate  case  is  not  formed  at  20  cm,  but  is  well-formed 
at  35  cm.  This  resu.l+-  should  bring  attention  to  the  fact  that 
ionizing  shock  strv=ii^Lii  plays  an  important  role  in  shock  for¬ 
mation  time  in  electromagnetic  shock  tube  technology. 
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SHOCK  FORMATION  STUDY 
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TIME  »  Zfis  /cm 


